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CN . Abstract 



The production and the final abundance of gravitino dark matter appear to depend crucially on the restoration 

^^ • of the global U{1) _R-symmetry of GMSB sectors in a threefold way. An _R-symmetric phase effectively suppresses 

^^ I the production of goldstinos from scatterings with the supersymmetric Standard Model particles and it generically 

^\ . initiates the goldstino production from the thermalized messenger particles. In addition, the GMSB spurion gets 

^N) ' displaced from the zero temperature minimum and under certain conditions it dominates the energy density of 

the universe producing late entropy. We show that it is possible to have high enough reheating temperatures that 

thermal leptogenesis and thermal vacuum selection can be realized without gravitino overproduction. The gravitino 

dark matter can be produced either thermally or non-thermally. In the former case the messenger scale has to be 

less than about 10^ GeV with the gravitino relatively heavy, 7713/2 > 0(10) GeV. In the later case, the gravitino is 

Ph' generically produced by the decay of the GMSB spurion field a process that always takes place for large messenger 

5s ' scales. A connection of our results with current collider and observational data is performed. 
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1 Introduction 

Supersymmetry is the most motivated and studied candidate theory beyond the Standard Model of particle 
physics. In supersymmetric models there is a fermionic degree of freedom for each bosonic one. This 
symmetry addresses the hierarchy problem between the electroweak scale and higher scales like those of the 
Grand Unified Theories (GUT) because it provides a natural mechanism for canceling divergences. Striking 
predictions of supersymmetric theories are the gauge coupling unification and the presence of particles with 
weak or sub- weak interactions that can account for the dark matter of the universe. Experimentally, we 
have not yet found any superpartners of the observed particles or profound supersymmetric deviations 
from the Standard Model. However, the discovery of a Higgs-like scalar at the Large Hadron Collider 
[1] [2] urges upon the hierarchy puzzle. The experimental data of the following years are expected to be of 
decisive importance for unrevealing the physics and symmetries of the TeV scale. 

If the physics beyond the Standard Model is supersymmetric it is necessary to understand the effects 
of finite temperature in supersymmetric theories because they are expected to describe physics in the 
early universe. The supersymmetry introduces scalar degrees of freedom hence the vacuum structure of 
the theory becomes complex. The thermal effects can either drive phase transitions in the MSSM or in 
the gauge mediated supersymmetry breaking sector [U IH [5l |6]. Also, if the early universe experiences 
temperatures larger than the sparticle masses then the supersymmetric degrees of freedom will get excited 
affecting the evolution of the universe and if longlived can be part of the dark matter. In supersymmetric 
theories the Lightest Supersymmetric Particle (LSP) can be stable. When the supersymmetry breaking is 
mediated by mainly Standard Model gauge interactions the LSP is generically predicted to be the gravitino 
and its present relic densisty has to be Q,^/2h^ ^ J7dm^^ — 0.118 [7^ where h is the Hubble constant in 
units of 100kms"^Mpc"^ 

The gauge mediation supersymmetry breaking (GMSB) scenario is well motivated from the particle 
phenomenology perspective. It is basically a calculable theory and provides a compelling explanation for 
the absence of excessive flavour-changing phenomena [5]. In addition to the phenomenology, it can be 
a cosmoUogically viable theory. It is free from Polonyi-like and gravitino problems due to the charged 
messenger sector that lies in an intermediate scale Mew ^ Afmess ^ Mp\. Actually, in the GMSB context 
the gravitino is a "problem" in the sense that the dark matter energy density observed in the universe can 
be explained in terms of thermally produced gravitinos only if a tuning is applied either in the reheating 
temperature (for not-equilibrated gravitinos) or the gravitino mass value (for gravitinos equilibrated with 
the MSSM fields). In the later case the gravitino mass has to be ~keV range and this value is disfavoured 
by the cosmological data [71 [9] . 

The gravitino is by definition a totally hidden sector sparticle. If the gravitino had only gravitational 
interactions the gravitino cosmology would be simple. However, this is true only for the helicity ±3/2 
component. Once supersymmetry breaks the gravitino becomes massive via the supcrhiggs mechanism and 
a goldstino fermion accounts for the ±1/2 helicity component. For gravitino LSP the goldstino couplings 
with the MSSM are enhanced compared to the ±3/2 helicity component. Moreover, the goldstino fermion 
can generically be coupled with other hidden sector fields. These imply that the gravitino dynamics may be 
subtle and the calculation of the final yield non-trivial. In most of the cases it is the reheating temperature 
or a temperature related to the hidden sector mass scales that controls the final yield albeit, it is very 
possible that non-thermal processes can be the dominant source. Apparently, the gravitino LSP cosmology 
is radically different with the neutralino case which accounts for a visible sector sparticle that is part of 
the thermal equilibrium until the, comparetively low, freeze out temperature Tq = 0(10 GeV). 



The calculation of the gravitino final yield is a compound task that involves knowledge of the MSSM 
properties, the behaviour of supersymmetry at finite temperature and the secluded/hidden sector dynamics. 

1.1 Overview 

The following superpotential is the basic paradigm of gauge mediation supersymmetry breaking: 

W ^FX + XX(j)4> , (1) 

where X is the Standard Model singlet spurion that breaks the supersymmetry, 0, the Standard Model 
charged messenger fields that mediate the supersymmetry breakdown to the observale sector and A is the 
messenger coupling; for review see ref. [8]. Here we are working in the spurion limit, where the effects of 
supersymmetry breaking are encoded by the expectation values (X) = Xq + O^Fx and any dynamics of the 
X field can be neglected. Generally, we also assume the messenger sector is weakly coupled. This model 
has a particular global U(l) symmetry, the i?-symmetry, under which R[ip] — R[(l>] =0 and R[X] = 2. 
This symmetry is important in models of spontaneous A'^ = 1 supersymmery breaking |10[ 111) . It is an 
accidental symmetry, since it has not been imposed on the theory but it is consequence of supersymmetry, 
gauge symmetry and field dimensionality. In order for the superpotential term J d'^x[W]p to conserve R, 
the superpotential itself must have R[W] = 2. The i?-symmetry is broken spontaneously by the lowest 
component vev of X. However, the Xq is not determined at the tree level. In the case of canonical Kiihler 
and in the tree level approximation there are only supersymmetric vacua and the X is a pseudomodulus, 
i.e. it is a flat direction with constant potential value F'^. The Coleman- Weinberg potential from the 
interaction between messengers and the spurion usually results in minimum at the origin that is not stable 
at the messenger direction. The spurion can be stabilized due to corrections to the Kahler potential for 
the spurion. A general expression is of the form 

K=\X\' + eJjl + eJ-^. (2) 

where A* is a cut-off scale of the theory related with the microscopic structure of the theory. For 64 = 1 
and eg < the minimum lies at Xq ~ A»/-\/|e6| and the i?-symmetry is spontaneously broken [T4j [151 E]- 
In the case that there is a bare mass in the messenger sector apart from the vev of X the i?-symmetry 
may break explicitly. For example, in the superpotential ([T]) we can add an i?-violating messenger mass 
5W = M(t)(j). A dimensionful constant, 5W = c, that cancels the potental energy of the phenomenologically 
acceptable vacuum is also an i?- violating contribution. For e4 = — 1 the minimum in the X direction lies 
at Xo = y3A2/(6Mpi) when 6W = c^ or, Aq = -M/A when 5W = Mcj)^ [13]. 

Majorana gaugino masses are closely related to the vacuum structure of the theory, see e.g. ref. [Tn[T6l 
I17j . If supersymmetry breaking is not accompanied by i?-breaking a split spectrum of gaugino and sfermion 
soft masses emerges. The i?-symmetry forbids the appearance of Majorana gaugino masses, m\, A- and 
/x-terms. The cosmological implications of the thermal restoration of the C/(l)ij symmetry were mentioned 
in ref. |18j . Namely, an i?-suppression of the mx and A- soft terms correspondingly suppresses the helicity 
±1/2 gravitino component production rate from MSSM thermal scatterings which is controlled by the ratio 
"^i/"^3/2- ^'^ i?-symmetric phase may also affect the gravitino mass. In supergravity a non zero value 
for the scalar component of the superpotential (W) 7^ is neccessary for tunning the vacuum energy to 
a nearly zero value and this explicit cancellation of the cosmological constant breaks the i?-symmetry, a 
breaking communicated to the visible sector by the gravitational interactions. It is this cancellation that 
gives a gravitino mass 

-3/. = e-/(-p>)i^, (3) 

where Mpi = 2.4 x lO'^^ GeV is the reduced Planck mass. However, generally in an i?-symmertic phase 
the ratio rn\/m^i2 is expected to be smaller than one. Indeed the scalar soft masses, m, contribute to the 
vacuum energy a fact that implies a lower bound for the gravitino mass while the Majorana gauginos mass 
can nearly vanish [T7] . 

The gravitino acquires mass by absorbing the massless goldstino fermion that appears once local super- 
symmetry breaks spontaneously. At high energies the goldstino corresponds to the helicity ±1/2 gravitino 
component. The GMSB goldstino has a tree level coupling with the messenger fields which can generate 



goldstinos more efficient than the MSSM fields. Indeed, for temperatures larger than the messenger mass 
scale, Minoss, the messenger fields number density have an unsuppressed thermal distribution and generate 
goldstinos via 2^2 scatterings with cross section scaling like A^/s, where s is the center of mass energy. 
This cross section peaks when T ~ M^^ss and is generally large enough to bring the goldstinos into thermal 
equilibrium regardless the gravitino mass value. Moreover, the messengers at temperatures T ~ M^^ss will 
decay into goldstinos with a width Fa oc X'^F'^/M^^^^. For sufficiently small messenger scales, Mmess ^10^ 
GeV, the goldstino yield from messenger decays (if this branching ratio is not negligible) can dominate 
over the one from scatterings. The conclusion is that once messengers acquire an equilibrium abundance 
the gravitinos get overproduced. Therefore, it seems that in the absence of any late entropy production 
the reheating temperature, T^h, has to be less than the messenger mass and the natural, thermal selection 
of the supersymmetry breaking vacuum [6] as well as the thermal leptogenesis scenario [19] cannot be 
realized. 

A small value for the superpotential coupling A does decrease the final gravitino yield from messenger 
scatterings and decays. However, in order the dark matter constraint J73/2^^ < 0.11 to be fulfilled the 
messenger scale has to be small and the coupling ultra small, A < 10~^°. Nevertheless, there is also another 
way to suppress the interaction strenght between the GMSB messenger fields and the goldstino fermion. 
The messengers couple with the X superfield that may not be the only source of supersymmetry breaking. 
It is a possible scenario that the goldstino does not reside in a single chiral superfield. This is the case when 
the supersymmetry breaking scale is divided into the secluded GMSB and a hidden sector that mediates 
the breaking to the visible one only via gravity. In the same fashion with the spurion of GMSB we can 
parametrize the supersymmetry breaking effects of this hidden sector, by the vev of another SM singlet 
spurion superfield Z, (Z) = Zo + O^Fz- Therefore we have the relation 

F^ = Fl+Fl, (4) 

where F is the fundamental scale of supersymmetry breaking and Fx is the scale of supersymmetry breaking 
felt by the messenger paticles. The XFx is the mass splitting inside the messenger supermultiplets and 
its size dictates the critical temperatures of the secluded GMSB sector [6|. A direct consequence of the 
supersymmetry breaking scale partition is that the goldstino fermion becomes a linear combination of the 
the fermionic components ipx and V'z of the superfields X and Z respectively. 

The ratio k = Fx / F signifies whether the supersymmetry breaking is stronger in the gauge mediation or 
the gravity sector. Notwithstanding, we assume that the gauge mediation contribution to the SM sparticle 
massses is the leading one, hence, the gravitino mass 

is smaller than the radiatively generated soft masses. Otherwise the attractive features of the gauge 
mediation as the calculability and the absence of the flavour problem are lost. Actually this is a working 
assumption here for the gravitino is considered the LSP and dark matter candidate. 

When fc '^ 1 the spurion X is the dominant source of supersymmetry breaking and the basic gauge 
mediation features remain intact. On the other hand, when fc <C 1 the extra dynamics of this hidden 
sector have to be taken into account!^. Since this sector breaks supersymmetry as well we consider that it 
is also accompanied by a U{1)r symmetry that is preserved in the vacuum state. This way the Z-hidden 
sector vacuum becomes an enhanced symmetry point and thus attractive during the cosmological evolution. 
Moreover, the gravity-mediated contributions to the sparticle spectum increase mainly the sfermion masses 
and not the i?-violating gaugino masses. Another important remark is that in the A; ^ 1 scenario the 
coupling of the gauge mediation messsenger fields with the goldstino field are suppressed because messengers 
couple only partially with the actual goldstino. Therefore, the goldstino production from the messenger 
fields is much less efficient and the gravitino overproduction problem from the thermalized messengers, as 
we will present, ameliorates or even gets solved for low messenger scale. 

The estimation of the gravitino yield is controlled by the product of the messenger coupling A and the 
supersymmetry breaking parameter k. The observed dark matter density indicates that the A • k takes 
small values. However, the A • fc is related with the other supersymmetry breaking sector parameters F 
and Mmcss by the expression 

Afc = M^K (6) 



^Except if k is that small due to e.g. indirect communication of the supersymmetry breaking to the messenger superfields. 



and cannot become arbitrary small. The value of A determines the MSSM soft masses m ~ (a/47r) A and it 
is tied to the electroweak scale -though constrained by the LHC data A > 10^ GeV. For fixed A a vanishing 
A • k implies either a small messenger scale or a heavy gravitino. The parameter space that we consider 
here is the messenger scale to be Mmcss > 10^ GeV and the gravitino mass m^/2 < 100 GeV. 

Nowadays the beyond the Standard Model particle physics has entered the LHC era. The discovery of 
a Higgs-like boson in the mass region around 125-126 GeV, with a statistical significance at the level of 
5-sigma, by the CMS and ATLAS teams have a significant impact on the supersymmetric theories. It is 
actually challenging to explain this Higgs boson mass in the context of minimal gauge mediation. In the 
minimal setup, MSSM plus GMSB sector described by ([T]), the nih ~ 126 GeV can be explained by large 
yl-terms, above 1 TeV, generated through renormalization group evolution from high messenger scales |20j . 
Otherwise, the MSSM scalars are too heavy to be withn the reach at the LHC and the tension between 
the Higgs mass and MSSM naturalness deteriorates. Another possible scenario is the presence of marginal 
superpotential interactions between MSSM and messenger superfields. This way large A-teims can be 
generated directly at low messenger scales, see e.g. ref. [2T|[22]. 

It is an interesting result that both cases, with high and low messenger scale, that generate large A- 
terms can be cosmologically viable even if messengers are thermalized. The case that the messenegr scale, 
XXq, is high corresponds to relatively large values for the coupling A and goldstinos get overproduced. 
However, the spurion X causes a late entropy production automatically and the dilution can have, under 
particular conditions, the correct magnitude [34]. In the case that messenger scale is low the goldstino 
production from the messenger fields can be small enough that the goldstinos are genarated mainly by 
the MSSM fields without any dilution to interfere at low temperatures. We note that in the later case 
the gravitino production from the MSSM takes place at the characteristic temperature of the i?- violating 
phase transition denoted by T^. 

Therefore high reheating temperatures are found not to be problematic in the GMSB scenarios. On 
the contrary, high reheating temperatures are welcome both for thermal leptogenesis and for the thermal 
selection of the supersymmetry breaking vacuum ^ . 

1.2 Outline of the results 

The gravitino dark matter theory has been developed in several notable works. In the context of the 
MSSM the results of ref. [23l [24j [25j [26l [27l established the relation between the gravitino abundance and 
the reheating temperature. In the context of the GMSB models the messenger sector changes essentially 
the final gravitino yield from the thermal plasma. In the ref. j28) it was shown that for temperatures 
larger than the messenger mass the messenger sector contribution to the gravitino yield is much larger 
than that of the MSSM because the messenger particles have stronger couplings to the goldstino and thus 
a stringent upper bound on the reheating temperature was found -apart from particular choices for the 
masses 7713/2 and M,„ess- Moreover, for viable cosmology the messenger fields cannot be stable unless the 
lightest messenger particle has mass less than 5 TeV [29]. Hence the messengers, being unstable particles, 
can have a further impact on the cosmological evolution. If the lightest messenger decays slowly enough 
then it can be a source of a late entropy production diluting the gravitino and any pre-existing abundance 
[30] . The decay width of the messenger particles depends on the messenger number violating interactions 
and hence the lightest messenger dilution magnitude can vary significantly |31j . Another basic degree of 
freedom of the GMSB models is the spurion field. In the case that the spurion oscillations dominate the 
energy density of the universe then it will produce a late entropy among with non-thermally produced 
gravitinos. It is found in ref. [311 1331 131] that the yield of these non-thermal gravitinos is possible to be of 
the right order of magnitude. 

In this work we re-investigate the physics of the gravitino production in the thermalized plasma of 
the early universe. We take into consideration that thermalized messengers induce a thermal mass on 
the spurion and for temperatures larger than the T^, which is the characteristic temperature of the R- 
symmetry breaking, the spurion vev is shifted to the origin of the field space. Following the result of ref. 
[T8] we assume that for T > Tf^ the gravitino production from the MSSM plasma effectively vanishes. On 
the other hand the relativistic messengers have a thermal distribution and efficiently generate goldstinos. 
We observe, actually, that for not too small superpotential coupling value. A, the goldstinos generated by 
messenger scatterings acquire a thermal distribution. 

Motivated by recent phenomenological studies of GMSB models that suggest a low messenger scale with 
marginal superpotential couplings of messengers with the MSSM sector for interpreting the LHC data, e.g. 



ref. [21] [22], we focus on the low messenger scale scenario and we examine its cosmological implications. 
We find that for low messenger scale the gravitino can have an acceptable relic density, ^3/2 < ^dm, only 
if the production of goldstinos from the thermalized messenger particles gets suppressed. This happens 
if the gravitino is relatively heavy TO3 /2 > 1 GeV and the superpotential coupling is small enough. The 
latter finding is also desirable for the thermal selection of the supersymmetry breaking vacuum [B]. The 
former, i.e. a heavy gravitino, indicates that there is an extra sector that breaks supersymmetry, but 
not the C/(1)_R, and is not coupled with SM charged messengers fields. Thus the gravitino mass is k~^ 
times enhanced despite the low value of the Mmess- Moreover, the marginal superpotential couplings of 
messengers with the MSSM render the lightest messenger shortlived and hence its decay does not dilute 
the thermal plasma. Therefore, if the messenger scale is low then the gravitinos get generated mainly by 
the MSSM plasma at the temperature Tf^. The cosmological constraint for the gravitinos, ri3/2 < ^dm, 
can be satisfied in a part of the parameter space where 10'' GeV < Mmess ^ 10^ GeV, m3/2 > 1 GeV and 
reheating temperatures Tfi < Trh < T±^/2. The T'-|-3/2 stands for the temperature that the helicity ±3/2 
gravitino component, generated via gravitational interactions, saturates the f^DM bound. We note that 
values for the messenger mass M^^ss < 10^ GeV can be reached only for A < 10^ GeV which is a rather 
constrained possibility according to the LHC data. 

There is a variance with the results of ref. [28] which concerns the goldstino yield from the messenger 
decays. Here, we find that the branching ratio of the messengers to goldstinos is generically too small to 
have a leading contribution to 5^3/2 because we consider superpotential couplings of the messengers with 
MSSM superfields. 

Above, we have implicitly assumed that the spurion X will reach the zero temperature minimum 
without large oscillations. However, this cannot be always the case. Actually, we find that this generically 
happens at low messenger scales. For larger values of the messenger scale, Mmoss ^ 10^° GeV, the spurion 
oscillations always dominate the energy density of the universe and dilution takes place. We mention that 
in the scenarios where the lightest messenger decays slowly enough to dilute the plasma, as described in 
ref. [3D], the spurion afterwards reheats the universe and therefore, the messenger dilution effects can be 
generally neglected. 

There are also phenomenological GMSB models that suggest a large messenger scale in order to explain 
the 126 GeV Higgs mass [20]. These models predict large mass for the gluinos rrig > 3 TeV. However, as we 
explained, large masses for the gluinos do not deteriorate the gravitino cosmology due to the fact that in 
the large messenger scale case the gravitinos are produced mainly non-thermally by the spurion decay. In 
this work we point out the part of the parameter space where the non-thermally production of gravitinos 
takes place without carrying out any detailed calculation of the non-thermal gravitinos abundance. This 
can be found in the ref. [32j |33l |34] . 

The article is organized as follows. In section 2 we discuss the behaviour of the supersymmetry at 
finite temperature mentioning issues as the vacuum selection and the values of the soft parameters in the 
reheated early universe. We also mention the important effects of an i?-symmetric thermal phase reserving 
a more detailed discussion with explicit GMSB models for the appendix A. The section 3 contains the 
interaction terms of the gravitino and the goldstino superfields with the MSSM and the messenger sector. 
There we give the relevant cross sections and decay rates for the calculation of the final gravitino relic 
density. In section 4 the total gravitino yield from MSSM and messenger sector is computed taking into 
account possible late time production. This section contains the main results of this work. Complementary 
calculations and some technical details can be found in the appendix B. In section 5 we briefly comment 
on supersymmetry breaking theories where the i?-symmetry restoration can be possibly realized without 
the thermalization messenger fields; such theories would actually give a different gravitino cosmology. In 
section 6 we combine our results on the gravitino relic density with GMSB phenomenology of the LHC 
era. We also discuss complementary cosmological constraints as the large scale structure in the Universe, 
the Big Bang Nucleosynthesis (BBN) and the leptogenesis scenario. In the last section we conclude. 

2 Supersymmetry breaking at finite temperature 

At finite temperature the fields that interact with the thermal plasma are no longer in their vacuum state. 
The occupation numbers are given by the Bose-Einstein formula. The temperature dependent one-loop 
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Figure 1: The figure shows the thermal evolution of the spurion vev, X(T), over the zero temperature vev Xq for 
different values of the cut-off scale A* = 10~^, 10""^, 10^^ Mp\ (from top to bottom) and for superpotential coupling 
value A = 10~^. The model considered is the ([1]) plus SW — c and SK = \X\'^/Al. It is apparent that as the 
temperature increases the spurion thermal average value approaches the origin of the field space. 



effective potential is of the form [35l |36l |32l [5] 

V^ - T* dxx^\n^l± exp 



-^x^+Mf/T^ 



(7) 



where Mf is an eigenvalue of the mass squared matrix. 

Knowledge of the behaviour of the supersymmetric field theory at finite temperature is basic for reliable 
cosmological calculations. The gravitino cosmology is a representative case of a such calculation. Its relic 
abundnace depends on the precise values that the scale of supersymmetry breaking, the soft masses, the 
messenger scale and hidden sector couplings have in a thermal phase. A particular example is the thermal 
restoration of the i?-symmetry that appears to be a generic phenomenon in GMSB models. In the ref. 
[T5] it was shown that the gravitino production from the MSSM sector is fairly insensitive to the reheating 
temperature once i?-symmetry restoration takes place. In this work, we also include the contribution from 
the thermalized messenger sector to the final gravitino yield. Another important issue tightly connected 
with the gauge mediation models is the metastability of the supersymmetry breaking vacua. The vacuum 
structure of the supersymmetry breaking potential at finite temperature has been studied in [5] where the 
danger of transition to a supersymmetric vacuum was underlined. 

High reheating processes as the leptogenesis or thermal gravitino production with mass in the MeV- 
GeV range and ^^3/2 ^ 0.2 cannot be secure if the phenomenologically acceptable minimum is thermally 
unstable. 



2.1 Thermally safe models of supersymmetry breaking 

By construction the gauge mediation scenario addresses the Polonyi cosmological problem thanks to the 
Yukawa couplings of the spurion with the messengers; however, it is these particular couplings that render 
the supersymmetry breaking vacuum metastable . The existence of many competing vacua, supersymmetric 
and non supersymmetric, poses the question of how natural it is for the system of fields to settle down 
into the phenomenological relevant vacuum with broken supersymmetry. In ref. ;6^ it was shown that the 
critical condition is the supersymmetry breaking sector to be sufficiently weakly coupled to the messenger 
fields. In addition the spurion pseudofiat direction (at the tree level) has to be lifted by radiative corrections 
from interactions with fields other than the messengers; in other words, the messengers one-loop correction 
to the spurion tree level potential has to be a subdominant contribution. Generally the one-loop corrections 
can be included in the Kalher. In summary, for the messenger superpotential, W — \X(j)(j), the following 
vacuum selection conditions have been found [6]: 

1. Canonical Kahler for the spurion X (regardless the value of A): thermally "discarded" . 

2. X-dependent corrections in the Kahler and A ^ 1: thermally favourable. 
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Figure 2: The effective potential for the spurion field (left panel) and messenger fields (right panel) at the temperatures 
5Tfi,Tft,0.5Tfi,0.2Tg and T — from top to bottom. In the left panel, the breaking of the _R-symmetry is manifest 
for temperatures T < Tfi : X{T < Tfi) w Xq- At the same temperature, T^, in the right panel the second order phase 
transition has not yet taken place. We used the values A* = 10^'^A/pi for the cut-off, A = 10^^ for the couplind and 
F = IQ-'^^'GeV^ for the susy breaking scale squared. 

Thermally discarded means that, once the messengers get thermalized, the system of the fields evolves 
towards the supersymmetric minimum. In the thermally favourable case the supersymmetry breaking 
minimum is always more attractive. For the sake of comprehension we can demonstrate this interesting 
behaviour by assuming, as a gedunken experiment, that we heat a part of the universe up. If the heating 
temperature is less than the M^css then both theories, 1 and 2, behave in a similar manner. On the contrary, 
when the temperature exceeds the messengers mass the first class of theories will exhibit a thermal phase 
transition towards a supersymmetric state while the second class will exhibit a kind of "elastic" behaviour 
where the initial zero temperature state is restored after the cooling of the systenlj. It is thus illustrative 
to characterize the first class of theories fire-sensitive and the second class fireproof. 



2.2 Finite temperature effects on the supersymmetry breaking sectors 

Supersymmetry behaves differently from other symmetries. While most symmetries, with typical examples 
of gauge or global SU(N) symmetries, if broken at zero temperature they are restored at sufficiently high 
temperature, exact supersymmetry at zero temperature breaks at T ^ |38j . This is a consequence of 
different statistics for bosons and fermions [33] ■ If the auxiliary fields acquire a nonvanishing thermal 
average a zero-momentum fermionic (collective) excitation can be created with zero energy. This massless 
excitation couples to the supercurrent in analogy with the Goldstone theorem at T = 0. Since a Gold- 
stone fermion is associated with supersymmetry breaking at T 7^ 0, the breaking is spontaneous. The 
corresponding real-time fermion propagator at low temperatures has the dependence '40!] 



^(po,P = 0)~e-'"/^/7oPo. 



(8) 



This expression indicates that the Goldstone contribution vanishes as T — > for the residue of the Goldstone 
pole vanishes as T — ^ 0. At the one-loop level, at finite temperature, the F auxiliary field develops a 
nonvanishing thermal average, (F) a, that is proportinal to e~™' "^ x (T/m)"' for T -^ m, where m, is a mass 
scale of a massive Wess-Zumino model [40ll41j . 

In reference |42j it was shown that the thermal Goldstone boson does not mix with the zero temperature 
goldstino and hence, cannot play any role in the gravitino production. Furthermore, the authors of ref. 
[43] examined the leading thermal correction to the goldstino-matter coupling, go = TOq/|F|, and found 
corrections of the form mg(r) Ki ml{l -T'^\F\-'^ /2A) and \F{T)\ w |i^|(l - r2|F|-V48) where mo and F 
the zero temperature values. The conclusion is that the thermal corrections to the gravitino mass and its 
production rate are negligible and the zero temperature estimates are qualitative correct. However, these 
results are concrete only for temperatures T <^ yF. 

Due to the fact that the the gravitino yield from the MSSM sector appears to be dominated by the 
highest temperature of the thermalized plasma it is not suggestive to adopt an agnostic stance towards 



■^If we knew that somewhere in the observable universe collisions with (the ultra high) center of mass energy ^/s > Mmess had taken 
place and an instantaneous local thermal equilibrium could have been realized then we could discard the first class of theories for describing 
the physics beyond the Standard Model. Otherwise, a bubble with the supersymmetric global minimum could have been generated. 



the actual values of the masses and the couplings for temperatures above the ^fF scale. Here, in what 
concerns the soft masses at fininte temperature, we will consider the effects of the i?-symmetry restoration 
which takes place at T^, = AkF/{XA^^/N), a temperature not far away from the \/F scale. 

Goldstino - MSSM decoupling at high temperatures 

As we will expose in detail in the next section, the goldstino field, which at high energies accounts for 
the helicity ±1/2 gravitino component, couples stronger with MSSM gauge superfields rather than the 
chiral ones. In particular, the goldstino coupling to gauginos is proportional to m\T/F and dominates 
over the coupling to sfermions which is proportional to ra^J F. The thermal effects can induce extra and 
possibly dramatic modification to the goldstino coupling. First of all it is the suppression of the gaugino 
masses due to _R-symmetry restoration. Indeed, the finite temperature effective potential for the scalar 
component of the GMSB spurion superfield has a vacuum structure (illustrated in figures 1 and 2) that 
differs significantly from the zero temperature potential. It was found in ref. [TB] that for T > Tfi the 
minimum of the effective potential preserves the i?-symmetry a fact that implies vanishing gaugino masses 
according to the relation 

mxiT) X{T) _ 1 



mx X, ^^,_^ 



^ 6fl(T) . (9) 



Due to the zero temperature relation m\ = {a/4TT)Fx/XQ one may conclude that the gaugino masses, 
and generally the soft masses, increase like 1/X{T) as the X{T) — )■ . However, we expect that at finite 
temperature the soft masses do not to blow up as the X field approaches the origin thanks to the collective 
phenomena of the thermal bath. In particular, the soft masses are generated via the messenger fields 
that run in the loops. At the zero temperature the messengers mass is XX^ whereas for X{T) <^ Xq the 
messenger fields are relativistic and a plasmon mass ^ gT for messengers has to be considered. From 
another perspective, we expect that the interaction range of the relativistic messengers is not infinite but 
it is screened in a Debye fashion. Apparently, the exact value of the soft masses is a non-trivial issue and 
we will not enter into more details here for this is out of this work's scope. From these considerations we 
can infer that on the one hand, sfermion masses remain finite at high temperature due to the non-zero 
plasmon mass for messengers, - taking also into account a possible thermal corrected F-term for T > -yAF - 
and, on the other hand, the gaugino masses are suppressed due to the i?-synimetry restoration. Therefore, 
we can extend our results regarding the goldstino-MSSM couplings for temperatures higher than the T^ 
by neglecting the sfermion mass value altogether and considering that the gaugino masses nearly vanish. 

Furthermore, we consider the gravitino mass to remain roughly unsuppressed in an i?-synimetric phase 
because the superpotential value at finite temperature, W{T), does not vanish. These conclusions are an 
heuristic, though essential, rule for a consistent calculation of the final gravitino yield. 

Wc finally comment on the variance found in the literature concerning the goldstino-MSSM interaction 
for temperatures larger than the messenger mass. The authors of ref. [28l [34] argued that for energies 
T ^ M^oss one needs to consider diagrams in which the messenger fields are in external lines. At these en- 
ergy scales there is no direct local couplings of the goldstino to the MSSM particles at lagrangian level and 
the ipG decouples. On the contrary, the authors of ref. |31) claimed that in local supersymmetry there are 
tree level particle-sparticle-goldstino interaction terms, see eq. (TTlT) . hence the goldstino does not decouple 
from the MSSM for temperatures M^^ss S T and, thus, the estimation of the gravitino abimdance is sen- 
sitive to the reheating temperature. Here we follow the approach of the authors of ref. [31]. The effective 
lagrangian for the gravitino is obtained from the TV = 1 supergravity theory regardless the mechanism 
employed in the messenger sector to communicate the supersymmetry breaking to the visible sector |44] . 
The energies under considerations are much higher than the gravitino mass and the goldstino interactions 
can be readily obtained and, indeed, these include goldstino-MSSM contanct interactions. In the absence 
of any particular cancellation the goldstino abundance appears to depend on Trh- In our work, the new 
feature is that we take into account the i?-symmetry dynamics and, as we argued, the MSSM-goldstino 
couplings get suppressed above the T|j temperature. 

The messengers mass 

The behaviour of the messengers mass at finite temperature is crucial not only for the calculation of 
the gravitino yield from the MSSM sector but also from the thermalized messenger sector. In a ther- 
mal phase the messenger masses are temperature dependent, Afmcss = XX (T). It is interesting to notice 



that the messenger mass and the T^ do not coincide. There are two cases. Let us consider their hmits 
in order to emphasize on the differences: i) the T^ :^ Mmess and ii) Tfi <C Mmess- In the i)-case the 
messengers get thermahzed at the temperature T — Mmess and a thermal mass is induced to the spurion 
m\{T) = 4F^/Aj + X^T'^N /A which, however, is not enough to drive the spurion thermal vev to the origin 
of the field space. The temperature has to increase to the value Tft, in order to be X{T) ~ and R- 
symmetry restoration to take place. As the universe cools down first the _R-symnietry breaking takes place, 
the messengers become massive and afterwards, when T ~ Mmess, the messengers become non-relativistic 
and finally decouple from the thermal bath. In the ii)-case the thermalization of messengers simultane- 
ously brings the spurion vev to the origin. As the universe cools there is a subtlety: at the temperature 
T ~ Mmoss the messengers neither become non relativistic nor decouple from the plasma. On the contrary 
they remain massless. The temperature has to decrease more in order the messengers to become massive. 
This happens at the temperature T^ess — T^ ^ Mmess- 

When one calculates the final gravitino yield one should know whether the Tft is larger than the Mmess- 
Otherwise, if T^ <C Mmess, the final gravitino yield, as we will demonstrate in the following sections, may 
have been underestimated. 

A last comment on the potential behaviour of the gauge mediation theories above the messenger scale 
regards the microscopic structure. The gauge coupling unification in the GUT theories is an appealing 
feature of the supersymmetric theories and hence, the choice of the messenger representation is such that 
the couplings remain perturbative all the way to Mqut- Also we work in the spurion limit, where the 
effects of supersymmetry breaking are encoded by the expectation vaulues (X) = X + d^Fx- Nevertheless, 
these assumptions may actually not be valid. The spurion or/and the messenger fields may be compound 
states at low energies due to some hidden strong dynamics. In such a case, new physics appear above 
a certain scale which can essentially modify the goldstino production in the thermalized universe. For 
the sake of consistency, in this work we assume that all the gauge and superpotential couplings remain 
perturbative up to the reheating temperature, T^-h- 

3 Gravitino interactions 

The gravitinos are spin 3/2 particles and described by the Lagrangian for vector spinor fields. Its interac- 
tions can be read by the general N=l supergravity Lagrangian (see e.g. |45 [ l46 j ) which contains interaction 
terms between the gravitino field ?/)p and the Noether current of local supersymmetry J^ 

£^^ = ee^'^^-^^a^PpV'. - ee^/^ [W^*Vm'^''>- + W^^^.^r'^"^.] - (^^^f^J'" + ^-C-) • (10) 

The symbol e stands for the determinant of the vierbein e — det e^ and accounts for the spacetime 
curvature. For a supergravity Lagrangian with a chiral [f^.f^) and gauge supermultiplets (A^JjA") the 
leading dimension- 5 interaction terms read 

. (11) 
where K{f,f*) is the Kahler potential, Kij* the Kalher metric, W{f) the superpotential and 2?^ the 

covariant derivatives in the presence of spacetime curvature and gauge boson fields. 

In exact supersymmetry the gravitino field is massless and possesses only helicity ±3/2 modes. Sponta- 
neously broken supersymmetry gives rise to helicity ±1/2 goldstino modes and the gravitino acquires mass 
by the super- Higgs mechanism. For energies much larger than the gravitino mass and for TO3/2 not larger 
than the effective mass splitting in the supermultiplets the helicity ±1/2 components of the gravitino can 
be treated as the true goldstino modes. This is implied by the equivalence theorem and the ±1/2 gravitino 
component can be written for ^/s 3> W3 /2 

3to3/2 

where if^a denotes the goldstino. This is the case we study: TeV scale mass splitting in the observable 
sector supermultiplets with the gravitino the lightest sparticle with 7713/2 < 100 GeV. 
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^/'-n/oi^^/^^G (12) 



The goldstino field is the fermionic component of a chiral superfield that brealcs supersymmetry. Here 
this role is played by the spurion X which obtains a non-vanishing vacuum expectation value {X) = 
Xq + 9^Fx- Still, the X field may not be the only field that breaks supersymmetry. Other hidden sector 
fields, that we collectively denote them with Z, may have Fz 7^ and the total i^-term breaking value in 
such a case is F = \/F'x + F^. Hence, the gravitino mass reads 7713/2 = Fx/ky/^MY>\ where k = FxjF . 
Correspondingly, the goldstino fermion is a linear combination of the the fermionic components t/^x and 
■i/'z of the superfields X and Z respectively: 

^ X ^ y 

ipG = -prV'x + -yi'z- (13) 

We are going to consider both fc ~ 1 and k <^ 1 scenarios. 

3.1 MSSM sector 

Substitution of eq. (fT2|) into the expression (fTT]) the dimension-5 operators of gravitinos are converted into 
dimension-6 operators involving goldstinos. Namely, in the four-component notation and in the flat space 
time the gravitino interaction Lagrangian (|lip becomes 

/:e//. goldstino = ^-^=-— [7^7^]V'gA'^F;, -I- -^ ^^^Vg/k/* + h.C. (14) 



m3/2SVGMpi ''" "^3/2 VSMpi 

where the derivatives in the interaction Lagrangian were replaced by the soft supersymmetry breaking 
masses. It is manifest from the above Lagrangian that the goldstino decouples from the theory in the 
supersymmetric limit. For energies -^/s > mx, m; single gravitinos can be created or annihilated dominantly 
through processes that include gauge superflelds components x,y,z with cross section 

1 V- . , ^ 2.45?7n? + 9.251^1 + 26glml 
^tot = 2 L '^-^.^(^ + y^z + i,o)^ 24W, M2^ (15) 

where 771^ (A — g^w^h) the supersymmetric Standard Model gaugino masses and rjx^y equals 1 (3/4) for 

an incident boson (fermion). We notice that for m\ — > and rrif y^ the cross section vanishes i.e. the 

goldstino appears to decouple from the MSSM without having supersymmetry restoration. This is the case 

of i?-symmetry restoration. 

Since the gravitino is the Lightest Supersymmetric Particle (LSP) the heavier sparticles will decay to a 

Standard Model held plus a gravitino. The decay width into gravitinos is nearly the same for both gauginos 

and sfermions 

1 ~ 5 

r,,, (,^^^^)^__^_^ (16) 

where i = X, f. 

Apart from the interactions of the goldstino field with the Supersymmetric Standard Model (fTTI) it is 
crucial to include all the other interactions that may alter the final gravitino yield. Thereby, investigating 
the entire goldstino couplings is of central importance. 

3.2 Messenger sector 

The supersymmetry breaking is communicated to the visible low energy sector via the messengers super- 
fields. When these fields are charged under the Standard Model we have the ordinary gauge mediation 
scenario; when the breaking is communicated via gravity the Mpi suppressed interactions induce soft masses 
of order 7(7,3/2 ■ We require that the GMSB messengers come in complete, vector-like SU{5) representations 
and the SM gauge couplings remain perturbative up to the GUT scale. Messengers acquire supersymmetry 
violating masses from the non vanishing vacuum expectation value of the auxiliary component of the X 
field via the superpotential coupling XXcfxj). Due to this coupling the Lagrangian contains the following 
interaction term between the fermion ipx and the messenger fields 

-^^x^^ = -A(^xX0'/' + ^xX00 + h.c.) (17) 
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where x^j Xa the fermionic components of the (j), (f> messengers. The fermionic component of X can be 
written in terms of the goldstino fermion as ipx — [Fx / F)tpG + •■• and the above interaction, eq. (fTTl) . is 
cast in 

5C-G^4> = -AA;(V'GX00 + '>pGXct><l> + h.c). (18) 

In the minimal case of 5 + 5 in the messenger sector, the messenger superpotential at the GUT scale 
has the form Wmoss — AX55. Bellow the GUT scale, SU(5) is broken and the superpotential is splitted 
into Wmcss = A3X(/)303 + \2X(f>24>2- The (pz and <j)^ transform as an SU(3) colour triplet and the onshell 
Lagrangian density that describes its interaction with the gluons G°, hypercharge bosons B^, gluinos A| 
and binos A^ in terms of the four component Dirac spinor xa and Majorana gauginos reads 

5Cv4,x = *X37^^mX3 + (^^03x)^(i?M'^3,L) + (i?^03,K)t(^M03,i?,) + 
V2ff3 



%/2: 



!5i 






y3n,RhPRX3 + 2/30IlA^PlX3 + h.c. (19) 



where D^ = d^ — ig^f^G^ — iy^giB^ the covariant derivative with j/3 = —1/3, y^ = 1/3 and, Plm. = 
(1 ± 75)/2. The (ji^^L is the left squark messenger that corresponds to the (/)3 scalar of the messenger 
superfield and 03, _r the right squark messenger that corresponds to the ^g. In order to have the complete 
interaction terms between the messengers fields and the Standard Model vector superfields we have also 
to include the 2 + 2 messenger doublets. Then, the scattering cross section for single goldstino production 
of messengers 0, 0* and plasmons for energies -^s ^ Mmcss is given by [28] 

V [a(Amess + A'^^ss -^ SmSSM + iZ-g) + ct( Ancss + ^MSSM ^ A[^oss + V-g)] - 6 ^\^k^ - (20) 

^ — ' s 

A,A\B 

where Amess — 4>,Xcl> are messenger field components, Bmssm = A",V^ are gauge superfield components 
and ^ = (gj +3(72+8g|)/47r. The above cross section has an inverse dependence on the center of mass energy 
squared and it maximizes for s ~ Afmoss which is the threshold between relativistic and non-relativistic, 
i.e. Boltzmann suppressed, messenger particles (for further analysis see appendix B). The messenger cross 
section ([20|) has a different behaviour than the cross section of the goldstino production from MSSM fields, 
eq. ()15p , which is independent of the center of mass energy ^/s and yields a goldstino abundance with a 
linear dependence on the temperature. 

The lightest messenger particles may decay late or even be stable and contribute to the dark matter 
abundance. Their final abundance is determined by the thermally averaged annihilation cross section which 
for temperatures T <ti Afmess is given by 

(aann«) ^ ^ H dv V^ {o ,^^v) C- ^ ^^' ' ^ (21) 

2v2 Jo 
where Xmoss = M-mcss/T . Using the expansion CTann^ = a + bv^ + 0{v^) the expression (|2ip develops the 
form ((Tanni') — a + Qb/x + 0{x~'^). Namely, the thermally averaged annihilation cross sections to MSSM 
particles and goldstinos are parametrized like 

{a{^r ^ MSSM)«) = — ^ (a - -^) , (a(00* ^ ^Pg^g)v) = -^ (a' - -^) (22) 

-'''^mcss V 2:^mcss / -'"moss \ 2;incss / 

with A~B~3x 10-^ and A' = l/327r, B' = 15/12877 [13131]. 

The decay width of messenger particles to goldstinos is dominated by the superpotential coupling (|18p 
and reads 

rA(X0^0^G) = ^A4A:4-ZL, (23) 

mess 

Further couplings of the messengers with the other sectors and in particular the observable sector are 
possible and have been proposed in the literatute with either cosmological or phenomenological orientation, 
see e.g. '31' or 1221 respectively. Superpotential interactions are cosmologically expected in order to prevent 
the lightest messenger from being stable. Furtermore, a Higgs at 126 GeV it is hard to be explained in the 
minimal versions of gauge mediation plus MSSM; introducing marginal superpotential interactions between 
the messengers and the Higgses can ameliorate or even solve the problem. The cosmological implications of 
these extra messenger interactions, regarding the calculation of the gravitino abundance, will be discussed 
in the 4.2.2 subsection. 
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3.3 Further hidden sector structure 

In this work we also investigate the scenario that the GMSB sector has a subleading contibution to the 
fundamental supersymmetry breaking scale, F, which is mainly sourced by the Z- hidden sector i.e. F ~ Fz- 
Respectively, the goldstino fermion is mainly composed by a Z-fermion. Generally, the gravitational 
couplings between the gravitino and the other hidden fields Zi can be derived from the gravitino mass 
term C = -e'^/'^ {il^^iCrt"'^l;i, + h.c.) where G ^ K + \n{W*W) and we have set Mpi = 1, see e.g. 05]. After 
the Taylor series expansion of the exponential we get 



i(e^/2G,)z,V'^a''''V',. (24) 



^ 2 

If the number density of these hidden sector particles in the early universe is large enough and the branching 
ratio to gravitinos is significant then the Z-hidden sector can be cosmologically problematic. The study of 
non-thermal production of gravitnos from other than the GMSB secluded sector is very model dependent 
and out of the scope of this work. Here we assume that the Z-sector is naturally equiped with an exact 
C/(l)ij symmetry [TU] which is important for justifying that the Zi scalars are not displaced from their zero 
temperature minimum in the post-infiationary era. Therefore the SM gauge singlet scalar, Z, which has no 
Yukawa couplings to messengers may not cause a Polonyi-like cosmological problem and not overproduce 
gravitinos via its decay. 

4 Gravitino relics from the thermalized plasma 

The gravitino number density 713/2 in the thermalized early universe evolves according to the Boltzmann 
equation [35] 



dn 



^^ f3i7n3/2 = (7sc+7dcc)|l-^|. (25) 

The 7sc = 7^^^^ + 7sc°'^'' is the gravitino production rate in scatterings with thermalized Standard Model 
sparticles and messengers and reads 






7sc = 0.14^ II + "^^;^^/^ 1 + ^>^'k'-fis) , (26) 



for six effective quark-squark flavour multiplets and mg(^) = mg(T)g'^ (fx) / g'^ (T) with /i about the elec- 
troweak scale. The /(s) is a dimensionlcss function of the center of mass energy squared that can be found 
in the appendix. The 7dcc accounts for the gravitino production rate from decays of particles in thermal 
equilibrium 

N WmSSM A^mcss 

_V^ cq^Tl^p _ V^ eg "Ij MSSM , V^ „cq "^' j^iness /r,y^ 

Idcc -2_^n- — yidec- 2^ rij T^^dcc + 2^^l Tg)^ doc l^'J 

where rrii/ (E) is the thermal average of the time dilation factor, ttij ~ m the soft masses and mi ~ M^^ss- 
We mention that the messenger particles share roughly the same mass for small messenger mass splittnigs, 
i.e. -M^pgg ^ XF. The decay rate of the Standard Model sparticles into gravitinos is given by the eq. p^ 
and that of the messengers by the (p3)) . 

4.1 Thermal production of gravitinos 

4.1.1 Supersymmetric Standard Model sector 

In a supersymmetric thermal bath, where mgoft = 0, the gravitino possesses only the helicity ±3/2 compo- 
nents and interacts with Mpi suppressed interactions. A thermal distribution for gravitationally interacting 
gravitinos can be achieved only for unrealistic temperatures T > tL^ ^ O(Mpi). In a softly broken su- 
persymmetric thermal bath the longitudinal mode of the gravitino interacts with the observable sector 
with strength enhanced by m\/m^i2. This could result in the thermalization of the gravitinos by the 
a -\- h -ir^ c -\- tj) scatterings in the MSSM plasma giving a freeze out temperature 

Til ""'^ - 2 X IQi^GeV f ^g^) V ITeVy 
3/2 VlOMcV/ \ rrig J ^ ' 
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Figure 3: The gravitino abundance with respect to the reheating temperature for masses m^/2 =0.1 (left panel) and 
1 GeV (right panel) considering only the scatterings from the MSSM sector. The red dashed line corresponds to the 
abundnance of the helicity ±3/2 gravitino component and the blue line to the ±1/2 component. The hidden sector 
parameters were chosen so as the gravitino abundance to be 513/2/1^ — 0.1 for reheating temperatures larger than Trh 
and smaller than 10^^m~L GeV^. The plateau in the panels is due to the i?-symmetry restoration i.e. due to the 
suppression of the helicity ±1/2 component. 



However, this is true only if T!;% < Tg^ which is usually not the case, unless the gravitinos are very 

light. For gravitino mass above the MeV range the Tg is smaller than the T°% and the gravitinos 

do not equilibrate with the MSSM plasma. The gravitino abundance in units of the critical density from 
scatterings with the MSSM plasmons for reheating temperatures Trh > Tf) is given by the expression |18) 



„MSSM (th), 2 



0.1 



/ SrhTfl, 

VlO^GcV, 



\ /GeV\ (7n^ 
) \m3/2j UTeV 



(29) 



where 6rh ^ 3/2 for T,.h > lOT^, /i ^ 100 GeV for the running of gluino mass nig from the electroweak 
scale to the i?-symmetry breaking temperature. For T^ ~ 4fcF/(\//VAA*) the abundance reads 

2 
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(30) 



We see that when three parameters A, A*, and k combine to the value AA,/fc ^ 10'^"'^ GeV then the 
gravitino can account for the dark matter of the universe. The dependence of the ^^3/2 
reheating temperature is depicted in figure 3. 



MSSM (th) 



on the 



4.1.2 Messenger sector 

Despite the fact that the _R-symmetry restoration essentially decouples the helicity ±1/2 gravitino mode 
from the MSSM thermal plasma it signals the thermalization of the messenger sector. Thermalized messen- 
ger fields can efficiently exchange energy with the goldstinos. Indeed, from the single goldstino production 
eq. (|20p we find that for 



^30 = 



:= (e^X^k^T^/s) = e^X^k^T > H{T) 



(31) 



the goldstinos can be also singly destroyed and hence, a thermal equilibrium for the the helicity ±1/2 
gravitinos with the thermalized messengers can be achieved. This is a case much similar to the axions 
that can obtain a thermal abundance due to the pion-axion conversion processes. In the above relation 
the H is the Hubble expansion scale H{T) — 0.33.g* T^/A/pi and g*{T) the effective number of degrees of 
freedom; assuming a pair of messenger multiplets in the 5 + 5 representation of the SU(5) and the MSSM 
in thermal equilibrium then it is g*{T > M^css) — 270. Plugging in numbers we find that the goldstinos 
acquire a thermal equilibrium distribution for coupling values 



Afc> 



T 



PI 
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(32) 
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where a coefficient of order 0(1) has been neglected. The Tmoss is the temperature that verifies the 
equation T = Afmoss(r). At the temperature Tmess the Boltzmann suppression of the messenger fields 
number density initiates which results, on the one hand, on the termination of the goldstino production 
from scatterings with the messengers and, on the other, on the exponential suppression of the induced 
thermal mass for the spurion X scalar field. The Tmoss does not coincide with the zero temperature 
messenger mass Mmess = AXq. Actually, it is always T,„oss < -^^moss- 

For the sake of clarity there are some comments in order here: Firstly, we underline that the messengers 
mass, defined as the curvature of the effective scalar potential at the messenger direction, is temperature 
dependent (as also explained in section 2). Hence, messengers fields become non relativistic when T ~ 
Mnioss(T) = \X(T) where X{T) < Xq; the explicit dependence of the spurion vev on the temperature can 
be found in the appendix A. Apparently, when it is Tmoss < Tfi it is a good approximation to say that 
Tnioss — -Mmoss bccause at the temperature T^ the spurion vev has approached the i?-violating Xq value, 
{X{Tj/i)) = Xo/2. On the contrary when Tj/i < Tmoss it means that Tmess < Afmess- We mention that this 
is a special case that does not lead to sufficient damping of the X scalar oscillation about the Xq vev and 
dilution takes place. This is due to the fact that the thermal mass for X disappears before X{T) reaches 
Xq. We will comment on this case later on (see subsection 4.2.4) when we will scan the full parameter 
space in order to estimate the final gravitino abundance. 

In order to simplify the analysis and formuli we will use the Mmess instead of the Tmess -except if 
otherwise stated- for the temperature that signals the transition of messengers to the non-relativistic 
regime. Actually, the final gravitino abundance is sensitive to initial temperatures when no dilution by the 
spurion takes place and this happens when the approximation Tmcss — -^moss is valid. 

Coming back to the abundance of the goldstinos we see that when the condition ([32]) is true the 
goldstinos number density is the equilibrium one. At temperatures T ~ Mmoss the goldstinos decouple and 
the asymptotic value of the yield, ^3/2.00: is the equilibrium value at freeze out |35j 

Y3/2 00 = Y'% (T ~ Mmess) = ^^ -^ = 0.278-^ (33) 

i/2,00 3/2 V mess; 27r4 g_^(T) g^(T) ^ ' 

where (/off = 2 x (3/4) and the yield ¥3/2 is defined by the ratio of the gravitino number density 71.3/2 
relative to the entropy density s. In terms of energy densities, the gravitino thermal relic density over the 
critical one reads 



3y^--<^)-«=(Sf)(il) 



(34) 

We stress that at T < T^ the goldstino interactions with the MSSM gauge supermultiplets are turned 
on and therefore the goldstinos may stay thermalized. If this is the case then the goldstino freeze out 
temperature with the MSSM degrees of freedom can be smaller than the Mmcss value. However, for 
gravitino masses above the ©(MeV) scale the freeze out temperature usually cannot be less than Mmcss- 

Unless the condition ((32)) is true the goldstinos do not equilibrate with the messenger fields. Never- 
theless, the messenger fields can generate a significant gravitino population. The abundance of decoupled 
goldstinos from scatterings with the thermalized messengers peaks at T ~ Mmcss f^nd was given in ref. |i"25] . 
We write the expression for the yield and the relic abundance in a form that manifests the dependence 
on the coupling A, the parameter fc and the messenger mass Mmcss in transparent way rather than the 
gaugino mass whose value is irrelevant in this kind of processes (irrelevant in the sense that there is no the 
rn^/?77,3 ,2 coupling strength coefficient here): 

and 

^me^ss(sc)^2 ^ 3.8 X 10^ Y^f^^^^'^ (^) . (36) 

For A ^ 1.5 X 10^ GeV the above relic density takes the simple form ViI^^^^^^'^'Ip? ~ 0.3 x lO^Afc. Equally, 

utilizing the relation Afc = (47r/Q!)?TiAMmcss/(v3-^^pi"T-3/2) the abundance p6p can be recast in the more 
familiar-looking form: 

..incss(sc) 7 2 ' vjeV 



3/2 \ m I 



)(t^)(^)'- <-' 
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These expressions give the goldstino abundance from the scatterings of the messengers in the plasma. Some 
technical details for the calculation of the ^372**'' can be found in the appendix B. 

4.2 Non-thermal production of gravitinos and dilution 

Apart from the gravitino emission in scatterings of thermalized supersymmetric particles there is also 
a second source that contributes to the final gravitino abundance: the decay of the plasma ingredients 
plus decoupled species that have potentially dominated the energy density of the universe. Hereafter, we 
are going to examine and present the gravitino production from the decays of the MSSM particles, the 
messengers and the spurion. 

4.2.1 MSSM decays 

The contribution to the r23/2 from the Standard Model sparticles decays does not depend on the reheating 
temperature. It peaks at temperatures T ^ fh where m, a soft mass. Utilizing the decay width into 
gravitinos ([T6| the relic density of gravitinos from MSSM decays is found to be j45] 



^--^->h^ c. 3.4 X 10-^ ^Z- -— -— . (38) 



4.2.2 Messenger decays 

There is no basic model building reason for the messenger particles of gauge mediation supersymmetry 
breaking models to be unstable. Actually, if the messenger fields have no direct couplings with the MSSM 
sector then the theory conserves a global messenger quantum number that can render the lightest messenger 
particle stable and therefore relic of the early universe [29] . The stable messenger can have an acceptable 
thermal relic population only under particular assumptions a fact that makes the stable messenger scenario 
not attractive. Even if the lightest messenger is unstable it may have an impact on the cosmological 
evolution by the production of late entropy. The decay of the heavier messenger particles can also generate 
an important number of non-thermal goldstinos. Actually, the decay channels of the messenger particles 
influence the final gravitino abundance. 

The decay width of messenger particles to goldstinos is dominated by the superpotential coupling (fT8| 
and takes the value Fa = fc^A^F^/(47rM^(,gj,). Also in this case the goldstino production from decays of 
the messengers peaks at temperatures of order of the Afmoss ■ 

We recall that the messenger sector consists of a TV pairs of chiral supermultiplets 4> + (j) and each pair 
describes a Dirac fermion with mass M = Mmess and two complex scalar particles with mass squared 
M'^±XF. Taking into account the SU(2) D-terms the mass of the scalars within the same isospin multiplet 
are further split. We denote the ligthest messenger as (J)lm- The scalar (pLM particle is lighter than its 
fermionic partner hence it cannot decay via the superpotential interaction ()18|) to a messenger fermion 
and a goldstino. If the messenger sector has a conserved quantum number corresponding to an accidental 
global symmetry the lightest messenger will be a stable scalar particle. According to the analysis of [29] . 
if the lightest messenger sits in representations of SU(5) it is either the neutral or the charged component 
of a weak doublet in the case of 5 + 5 or a weak singlet with a unit of electric charge in the case of 
10 + 10. The (pLM can be also an SU(3)x SU(2)x U(l) singlet in the case of 16 + 16 representations 
of SO(IO). For the first two cases the lightest messenger decouples from the thermal plasma when it is 
non-relativistic and the yield is roughly Y^j^^ ~ 3.7 x lO^^''A/0J^J^,/lO^GeV which gives the abundance 
"fa., = 2-8 X 10»yfc.,M,,„/GeV i.e. 

If the 4>LM is a Standard Model singlet it decouples when it is still relativistic with yield 0(10^'^) and 
its relic density is Vtff,^^,jh? ~ 7 x 10^^Mnioss/(10^ GeV) given that the reheating temperature was about 
or higher than the GUT scale where, actually, the thermal equilibrium is doubtful. However, there is 
a crucial point here that has to be taken into account: the supersymmetry preserving vacuum is the 
attractive minimum of the finite temperature effective potential [5] since the temperature corrections for 
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the SM singlet 4>lm are negligible and cannot compensate the tachyonic directions during the reheating 
phase. For this reason, we are not going to consider the Standard Model singlet (J)lm case altogether. 

The lightest messenger can be a dark matter component without overdosing the universe only under 
particular conditions hence, messenger number violating interactions that allow the lightest messenger to 
decay are necessary. An extended study of the messenger couplings in a cosmological context has been 
carried out in ref. |31j . Generally, superpotential couplings between the messenger and the matter sector 
are not problem-free. Although the SM gauge symmetries allow a large number of superpotential cou- 
plings there are constraints from the proton decay limits, from negative contributions to squared sfermions 
masses generated at one loop level and from bounds on FCNC processes. Nevertheless, extra marginal 
messenger interactions seem to be significant from a particle physics phenomenological point of view as 
well. Indeed, the difficulty of the gauge infractions in generating appropriate contributions to Higgs sector 
soft parameters, given that rrih ~ 126 GeV, suggests that a complete theory of supersymmetry breaking at 
the weak scale should include both gauge mediation and additional coupling to MSSM and in particular 
the Higgs sector. In ref. [22] a recent classification of all possible such interactions can be found. Here we 
give the messenger-matter couplings which were firstly introduced in the ref. [47] and have the form 

W D ylHD<p2e^ + y\HDQ^h , (40) 

where 4>2 and 03 are lepton and quark like messengers, Hd, e, Qi are respectively SM down- type Higgs, 
right-handed lepton and quark doublet and the j/* are Yukawa couplings with family index i. Tachyonic 
sfermion masses are avoided for ^ |yp < 10"'^ [45]. Flavour changing neutral currents lead to even stronger 
constraints. Assuming conservatively that yl ^ and zero Yukawa couplings with the second and third 
families then there is the limit on this Yukawa coupling y] < (Afmoss/lO^GeV) ^Tj. The superpotential 
couplings (j40]) induce a decay width 

= l^ . (41) 

Unless the Yukawa couplings y are extremely small the decay of the messengers takes place fast without 
entropy production i.e. before the domination of the energy density of the universe by the messenger fields. 
Actually, the decay rate (pi]) can be larger not only from the Hubble rate at the temperature T ~ M^^ss 
but also from the pair creation rate of the messenger particles themselves. Indeed, when Ty > Tint = (nav) 
where a ~ a^/T^ the messengers particles decay more rapidly than the rate produced by 2 -f-?- 2 scatterings 
in the thermal plasma. This simply means that for T > Mmcss the messenger particles might have an 
equilibrium number density due to the inverse decays of MSSM particles to messengers, which occur with 
inverse decay rate Fij, ~ Fy, than due to the two-to-two scatterings. It is worth to notice that if Tiy <g; Ty, 
i.e. T (or CP) invariance was violated, then it might have been possible to significantly suppress the 
messenger number density in the plasma and hence, the subsequent goldstino production. However, there 
is no any (at least profound) physical reason for such a suppresion to happen. The inverse decay rate of 
the messengers has the following temperature dependence [35] 

i 1,(7 J _ i , I (^^_/r)3/2,^M_/T T < M_ 

Nevertheless, another kind of interaction might exist and, unless the superpotential couplings (PH)) are 
present, can induce the decay of the lightest messenger. If this interaction is weak enough then it can 
cause a late entropy production. One can assume a messenger-matter mixing due to a correction in the 
superpotential 5W ~ {{W) /M|j)5m5_f where the subscripts M and F denote the multiplet that belongs 
to the messenger sector and MSSM matter respectively and (W) ~ TO3/2Mp[. It was proposed in ref. [3D] 
and it provides a decay channel of (J)lm into a Standard Model lepton and a gaugino with decay rate, F^, 
estimated to be 

J. _ 92 "^3/2 



167r M„ 



(42) 



We have assumed small mass splittings between the messenger particles, i.e. M^g^g ^ XF, and for order of 
magnitude estimations we use approximately the Mmcss instead of the M^^j^^ for the mass of the lightest 



^We note, however, that superpotential coupUngs of the form II40I I with sufficiently large Yukawa couplings could block the evolution of 
the system of fields towards the supcrsymmetric minimum. 
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messenger. 

Decay widths hierarchy and goldstino yield 

The amplitudes of the various decay rates of the messenger fields are essential for the estimation of the 
final gravitino abundance. As we will argue, the Ty, Tm and Fa values not only determine whether late 
entropy production takes place but, also, their ratio controls the goldstino yield from the messenger decays. 
Comparing the messenger-matter decay rates (|4T|) and (j42|) with the messengers-goldstino rate ([23]) one 
finds that for Afc ^ 1 it is Ty ^ Tx and Tm ^ T^, as illustrated in figure 4. In particular, taking the 
Yukawa coupling y — O.lAfnioss/lO'* GeV such that the flavour changing neutral current constraints are 
satisfied then 

^ -^ 2 X 10-- (^y (^—X ri2!^V (43) 

and 

lA . 1 4 f ^V f 12:^ V (44) 

In other words, for Afc <C 1, the messenger particles are much more probable to decay to MSSM particles 
rather than to goldstino fields. The total decay rate of the messenger fields is Ftot = ^\ + ^y-\-^rm therefore 
the goldstino yield from the messenger decays is 

^me..(dcc) ^ B,„N^^YZ.,{T ^ M^e..) (45) 

where ^3/2 = r^/rtot < 1 and N^^ the messenger degrees of freedom that decay to goldstino and a lighter 
messenger superpartner. The goldstino yield from the decays is dominated by temperatures T ^ Mmcss- 
Actually, the Ftot is much larger than the Hubble rate when messengers become non-relativistic: 



„ _ Ftot 



..„.*m3™j~^'<'<wl='^'^ '"' 



In the figure 4 the scaling of the three different decay rates is shown for different messenger mass scales. 
It is manifest that for Afc < 10~^ it is Fj, + F,„ ^ Fa and the Y^^ ^'^ , according the expression (P5)) is 
smaller than the goldstino yield from scatterings (|35|) . Taking into consideration the relation ((32)) . i.e. the 
fact that the goldstinos acquire a thermal abundnace for Afc > 10~^(Miness/10*GeV)^/^, we conclude that 
we always have 

^mcss(dcc) T^mcss(sc) ^ ^cq / .-x 

^3/2 ^ ^3/2 - ^3/2 y^' > 

regardless the messenger scale (see also appendix B). 

Due to the fact that we explicitly consider messenger-MSSM interactions our result is different than 
that of ref. [55] and [21] which find that the goldstino yield from the decays of the messenger fields with 
mass Mmoss ^ 10^ GeV dominates over the goldstino yield from scatterings. Therefore, the goldstinos are 
produced more efficiently from the scatterings rather than the decay of the thermalized messengers and we 
will consider the ^Q^n^** "'^ as a small correction to the final gravitino yield. The K,'^"^** "'^ is important 
only when the messengers dominate the energy density of the universe before their final decay. 

4.2.3 Dilution from lightest messenger decay 

Let us assume, here, that there are no direct messenger-matter couplings in the superpotential i.e. Fy = 0. 
However, in order the messenger sector not to be totally isolated from the observable, (this would imply 
that the lightest messenger particle is stable) we take, in this subsection Fj^ = and F„j 7^ 0. Vanishing 
direct messenger-matter direct mixings makes the lightest messenger (J)lm longlived. Unless (J)lm is a 
Standard Model singlet and F^ is weak enough it will decouple from the thermal plasma at the temperature 
Tf ~ M^^j^j20. Afterwards, the energy density of the (I)lm, P4>lm = '2Tr'^g*{T)T^MjncssY^^j,j4:5, wih scale 
like a~^ contrary to the a~^ scaling of the radiation, where a is the scale factor of the universe. The 
temporary relic messengers 4>lm become the largest component of the Universe energy density at the 
temperature 

^f^sF,,,, ^ 5 X 10-^ (-^^ ' GeV, (48) 
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Figure 4: The three decay rates of the messenger particles to goldstinos (red dashed hne -A) and to MSSM particles 
(continuous lines, blue -y for superpotential Yukawa couplings (UHl) and green -m for the correction (HU). The three 
panels correspond to different messenger masses, Mmoss = 10^ , 10^ , and 10^ GeV (from left to right). It is manifest 
that for fcA ^ 1 the messengers dominantly decay to MSSM degrees of freedom. 



where, for simplicity, we used the approximation M„ 
messenger-matter mixing [30] at the temperature 



Ma, 



The. 



>LM 



will decay due to the "indirect" 
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If T0J^J^J ^ Tdom then the thermal relic (^lm dominate the energy density, i.e. ptot — Prad+Pm ^ Prad where 
Prad — 7r^ff*(T)T'*/30, and its subsequent decay will produce significant entropy, reheating the Universe 
and diluting any pre-existing abundances. The entropy release is given by the dilution factor Amoss 



A, 



Ta, 



dom 



= 4.6 X kVx 



\ 3/2 

108 GeV ) 



A/„ 



(50) 



By definition, the dilution factor cannot be less than one. Apparently, there is no dilution from the lightest 
messenger decay if the messenger mass, the messenger superpotential coupling A and the hidden sector 
parameter k are small enough. Replacing the 7713/2 from the relation ([B]) we find that there is no dilution 
from messenger decays, i.e. Tdom < T'0J^J^^ , when 



Afc2 < 38 



lO^GeVy / A \ 



Mmcs J [lO^GeVj 



(51) 



Figures 5-7 illustrate this condition for different messenger scales. The absence of dilution may not be 
problematic because, firstly, the goldstino abundance from scatterings in the plasma, given by the expression 
(P5|) . has a dependence Kf^^'^''^'"^ ex A^fc^ and it can be small enough. Secondly, the spurion may be also a 

source of late entropy production. For \k > 10~^(Mnioss/108GeV)^/^ the goldstinos get thermalized and 
after the late entropy production the gravitino relic density obtains the value 



"3/2 

^mcss 



0.12 



1 



108 GeV y^' 

A3/2fc2 V Mmcss J 



(52) 



The dilution of the gravitinos from the lightest messenger particle decay induced by the superpotential 
coupling SW ~ {{W) /Mp)5m5f was proposed in ref. [30]. However, there is a small difference: here the 
goldstinos acquire a thermal abundance due to the scatterings of the messenger fields for any temperature 
T > Mmcss regardless the gravitino mass value, whereas in ref. [301 tbe goldstinos thermalize due to MSSM 
fields, i.e. the messenger contribution was ignored. 

4.2.4 Spurion Decay 

A basic advantage of gauge mediation schemes is the absence of the Polonyi-like problem. The super- 
potential coupling SW = XX(j)(f> decreases the lifetime of the spurion rendering it harmless for the Big 
Bang Nucleosynthesis (BBN) processes. Nevertheless, the spurion decays slowly and if there is a scalar X 
condensate that dominates the energy density of the unverse then dilution and entropy production take 
place. 
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The i?-syninietry restoration happens when the vev of the spuion, X, becomes vanishing. Thus in an 
i?-symmetric thermal phase the spurion finds itself dispaced from its zero temperature vacuum expectation 
vaule Xq. When the temperature falls to the Tfi value the minimum of the effective potential moves to Xq 
and the scalar field X follows the minimum without sizeable oscillations only if the thermal mass for the 
X field is not too weak [BJ [31] . Let us be more specific here. Taking the second derivative in the spurion 
direction of the finite temperature effective potential, the effective mass for the spurion was found to be [5] 

m3,(r)^yJ^^4^ + 7V^. (53) 

At the temperature T^ — 4F/(AA,) the thermal correction to the spurion mass squared, NX^T"^ /A, becomes 
equal to the zero temperature mass squared. The spurion field follows the minimum which moves, roughly 
instantaneously, towards the zero temperature vacuum expectation value when its mass is larger than the 
Hubble scale. At the temperature T^ the condition mx > H translates into 

Taking the ([T]) model with 5W = c as a representative example for an order of magnitude estimation of 
the above constraint on the coupling we find 

^S'°"(to4v)"'(t54v)"^' (^=) 

Although the constraint (|54)) is necessary it is not sufficient. One should take into account that the 
thermal mass for the spurion vanishes altogether when the messengers become non-relativistic (the X-field 
"decouples"). Hence, it has to be Tf- > M^essiTf-) = AX(T^) = \hR{T^)Xo = M,ness/2 (see appendix A) 
which translates into the bound 

.,...-.(i£-X)-(^)-. 

Unless constraints (IM|) and (I55|) are both satisfied then the energy stored in the spurion oscillations dom- 
inate the energy density of the universe and its subsequent decay will dilute the pre-existing abundances, 
see figures 5-7. The dilution magnitude, Agpur was calculated in the ref. |34| and was found to vary from 
^spur = 1 to Aspui- = 10^. The gravitino yield generated by the spurion decay, i^oT,""^, defined by the ratio 
713/2/5 where s the entropy density produced can be evaluated as 

•c^spur _ "'3/2 __ O pspur Jspur , .„■, 

^3/2 — ^^ - o-°3/2 "Z ■ \^') 



'3/2 - g 2^/2 



mx 



The Tgpur stands for the spurion decay temperature and B^J'^ for the branching ratio of the two gravitino 
decay mode. They are calculated in ref. [32j|33] to be Tspur = 0(0.1-10) GeV and Bf^l' = 0(10-10-10-^). 
Summing up, the gravitino relic density after the reheating of the universe from the spurion decay reads 

n,„ = -^+m,„Y^J^^^^ (58) 

^spur ' /Ocr,0 

where ^.f^^ = f}"^^ for Afc > (M^ess/Mpi)^/' or ^f^^ = 0^/7^'''^ + ^^s^m ^ ^grav otherwise. The 
pcT,o/so ~ 3.6 X lO^^h^^ GeV is the critical density divided by the entropy density at present. There are 
parameter regions where the dark matter is mostly non-thermally or a mixture of thermally (diluted) and 
non-thermally produced gravitinos j34] . 

4.3 Total abundance of the thermally and non-thermally produced gravitinos 

Here we sum up the gravitino relic abundance from the MSSM and the messenger sector. The final result 
depends on the reheating temperature and values of the coupling A and the hidden sector parameter k. 
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A) Thermally produced helicity ±3/2 gravitinos 

The helicity ±3/2 component interacts gravitationally with the thermal plasma. It can be produced 
mainly from 2 — > 2 scatterings with thermalized messenger and MSSM particles. The cross section is 
independent of the temperature (apart from the running of the gauge couplings) and hence the relic 
abundance has a linear dependence on the temperature. This contribution is always present: 

B) Thermally plus non-thermally produced helicity ±1/2 gravitinos from MSSM sector 

1. Trh > maxlMmcss,?"!?} 

^MSSM (sc)^2 ^ 15 X 1^ /10^:^GeV\ /m^^ 

^^3/2 ^ ^■^^'' ^\ \Klk )\lTc\) ^^^> 

2. m/10 < r,h < max{M„,ess, T^} 



3. Trh > m/10 

^MSSM(.ec),. ^ 3 , ^ ^,_. /Gc^\ / J28 y- /^y 



3/2 / - ^ 3 



C) Thermally plus non-thermally produced helicity ±1/2 gravitinos from messengers 
i.e. Trh > Mmcss 

1. Thermal goldstinos, Afc > {Mmcss/MpiY^^: 

VGcV/ \g*{T) 

2. Thermally produced goldstinos (from scatterings), Afc < (Mmoss/Mpi)^'^ 



n- .^ ~ 0(5) X 10- U-^ -^ (63) 



0»r'"'.' ~ 3.A^P (^) (^2^) (-|;L) ,«,) 

Summing up 

• Total gravitino abundance (when spurion's oscillations are adiabatically suppressed): 

n,/2{T,^) = min {nf;r + ^W^'"^ + ^f/7' ^1%} (es) 

when no dilution from the lightest messenger decay takes place. The gravitino abundance is diluted by 
the lightest messenger decay, i.e. ^^3/2 — > ^3/2/^mossj only under the special conditions: Fj, = and 
XP > 38 (lO^GeV/Mmcss)^ (A/lO^GeV). 

D) Non thermally produced goldstinos from the spurion decay and spurion dilution 

{ Trh > maxjAfmess, 7i?| 1 . .,, . , , 

^^ ^m/,^ 1, r <^m f ^ spuriou oscillatious arc undamped 

\A<5Tu/Mpi or M^ess > 2T^ J 

Hence, 

• Total gravitino abundance (when spurion's oscillations dominate the energy density): 
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Figure 5: Contour plots of the helicity ±1/2 gravitino abundance for different values of the superpotential coupling 
A, the parameter k and low messenger mass Mmoss- The gray region (lower left corner) corresponds to gravitino mass 
'713/2 > 100 GeV where gravity effects play a dominant role. In the rest contour area it is TO3/2 < 100 GeV. The light 
coloured region above the green line (upper right corner) corresponds to gravitinos with thermal abundance. In the 
white part of the contour the gravitino is not thermalized, albeit it is overabundant due to messengers scatterings. The 
yellow strip corresponds to A • fc values that render the gravitino production from messenger fields inefficient thus, the 
gravitino production is dominated by the MSSM gauge superfields; the gravitino relic density W3/2 = $13/2/1^ = 0.11 
can be realized. In the two area between the upper corners of the axes and the dashed lines dilution by the spurion 
decay takes place. 



f^3/2(Tl-h) 



{Qmess 
"3/2 



qMSSM 
"3/2 



Qgrav 

"3/2 ' 



^spur 



"3/21 



^ V {GeYJ 



(66) 



If there is Amoss > 1 then the final dilution factor is determined solely by the Agpurion since the spurion de- 
cays slower than the lightest messenger (though the T0^J^J is more model dependent). Actually, it happens 



that the parameter regions that implement Amoss > 1 and A, 



spurion 



> 1 largely overlap (figure 7). Thus, 



except for special cases, the potential dilution effects of the lightest messenger decay can be neglected. 



The total gravitino abundance is illustrated in the contour plots of figures 5, 6 and 7. All the plots give 
the 03/2/1^ for reheating temperatures higher than the messenger mass. Otherwise there are no thermal 
excitations of the messenger fields, i?-symmetry is not restored and the goldstino production is sourced 
only by the MSSM gauge superfields at Trh, i-e. the rehc abundance is given by the expression (pT|) . 

Generally, a thermal messenger population overproduces goldstinos apart from a parameter area (with 
yellow colour) where it is Afc ^ 1 and the messenger scale is low, Mmess < 10^ GeV. In the yellow areas 
of panels (a) and (b), figure 5, the goldstino yield from messenger scatterings is inadequeate to explain 
the observed relic density of the dark matter. Hence the messenger contribution to ^^3/2 is subdominant 
and the gravitinos are mainly generated by the MSSM superfields. The reheating temperature can be very 
high -the only constraint comes from the helicity ±3/2 component- and the gravitinos are generated at 
the moment that the i?-symmetry breaks. The gravitino relic density is given by the expression (j60p . The 
parameters A, k and A^ can be chosen such that fl^^h'^ =0.11 even when the gluino mass is rrig > 1 
TeV. Due to the smallness of the A • k product the gravitino mass, TO3/2 ^ 10~-^''M„iess/(Afc), is large: 
1 GeV < m,3/2 ^ 100 GeV. We mention that gravitino mass values 7713/2 < 10 GeV can be reached for 
messenger mass scales Mmess < 10^ GeV which are rather constrained by the LHC data. 

In the greatest part of the parameter space of the contour plots (a) to (f) the gravitino is overabundant 
(mostly thermal) due to the messenger scatterings. However, this parameter space is potentially not fully 
excluded. The entropy production due to the spurion decay makes a large part of the parameter space 
viable. There, the dilution by the spurion takes place naturally since it is offset because of the i?-symmetry 
restoration and its oscillations are not damped. The spurion oscillations are sizable at the left side of the 
red dashed line due to smallness of the spurion mass, mx{Tfi) < H(Tfi), and at the right side of the 
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Figure 6: The same with figure 5 contour plots for intermediate messenger mass Mmess- There is no (yelfow) area in 
this contour where ^^3/2/1^ < 0.11 is possible without late entropy production . The overabundant gravitinos can de 
diluted by the spurion decay which takes place in the two areas between the upper corners of the axes and the dashed 
lines. The correct magnitude of the dilution takes place for a specific choice of the parameters. 
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Figure 7: The same with figure 5 contour plots for large messenger mass Mmoss- Also in this contour there is no 
(yellow) area where n^/2h^ ^ 0.11 is possible without late entropy production. Here, there is another possible source 
of dilution: the messenger fields apart from the spurion. The yellow line in the upper right corner gives ^^3/2/1^ = 0.11. 
However, dilution by messenger decays requires special couplings and it overlaps the area where generically spurion 
dilutes. Due to the large messenger masses, we mention here that (all, a-f) contours give the helicity ±1/2 gravitino 
abundance for reheating temperatures T^\^ > Mmoss- H T'rh < Mmoss the gravitino abundance is simply given by the 
expression (pT|) . 
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blue dotdashed line due to the disappearance of the non-relativistic messengers from the thermal plasma, 
Mracss{Tft) > Tfi. A tuning of the parameters is required in order the gravitino abundance, given by the 
expression (|66|) . to be the correct one. According to [32l |33j |34] there are regions in the parameter space 
for gravitino mass between 10 MeV and 1 GeV where the observed dark matter density can be explained 
by the decay of the spurion. It is also possible that a dilution caused by the decay of the lightest messenger 
to take place albeit, as we commented, this is a rather restricted case. 

5 Separate mediation and /^-breaking sectors 

The models that we have studied until now share a distinct feature: the i?-syinmetry restoration takes 
place only if the messenger fields have thermalized. However, it is possible the sector which mediates the 
supersymmetry breaking not to coincide with the sector responsible for the i?-symmetry breaking. This is 
the case if the messengers masses are not controlled by the vev of the spurion field. In particular for X 
independent messenger masses and Mmess ^ Tf^ the goldstino abundance from thermal processes may be 
given solely by the MSSM sector. We can sketch an example of a model which implements such a scenario. 
The superpotential 

W^FX + \X(t>,<P2 + AX0201 + Mmcss (-^101 + 02<^2) + SWh{X, ...). (67) 

breaks supersymmetry due to the nonzero vev of Fx component of the spurion field, it mediates the 
breaking to the observable sector via the (j)i and (f>i fields which are charged under the SM gauge group and 
finally breaks the i?-symmetry via the non-zero vev of the spurion i.e. (X) = Xq + 9^F with XXq <C M^^^ss- 
The SWfi includes the hidden sector dynamics that stabilize the spurion in the Xq ^ _R-violating vev. At 
finite temperature the spurion field can be driven to the origin, X{T > T^) = 0, due to the 6Wh dynamics 
and not due to the messenger fields that stay out of the thermal equilibrium for Tih < Mmoss- If the 
goldstino production from the possible hidden sector thermalization is suppressed then the goldstinos are 
generated at the temperature Tfi by the MSSM vector superfields with relic density given by the expression 
([29|) . see figure 3. However, we should mention that the demo model ([67|) generates a "little" hierarchy 
between the sfermion and gauginos soft masses due to the fact that not all the messenger fields couple with 
the spurion field X. 

6 Collider and cosmological constraints on the gauge mediation 
parameters 

6.1 LHC 

The touchstone of the TeV scale supersymmetry is the ground based experiments. The Large Hadron Col- 
lider is nowadays the leading experiment shedding "light" on distances less than 10""'^^ meters. Positive or 
negative experimental evidences for spartiles can support or constrain further the supersymmetric theories. 
In the framework of specific models these experimental results can, in turn, be translated into constraints 
on the mass and interactions of gravitino which is an undetectable particle - thus a dark matter candidate. 
The LHC can indeed discover or give hints for the nature of the dark matter. The dark matter might be 
produced in the collider and, being stable and (sub)weakly interacting, it will carry away a substantial 
amount of energy, the so-called missing energy. The gravitino dark matter, being the LSP, can be pro- 
duced by the NLSP decays. Measurements of the NLSP lifetime and mass can actually reconstruct the 
gravitino mass, 7713/2, and the supersymmetry breaking scale X]j(l (^i) P)^ ■ According to eq. (fT6|) the 
NLSP lifetime is given roughly by 

CTNLSP ~ 487r— r^ = 167r-r-g 

"^NLSP "^NLSP 

Collider-stable NLSPs imply a relatively heavy gravitino while prompt decaying NLSPs a light one (ap- 
pearing as missing energy). 
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Despite the entire absence of sparticle signals some assumptions about the messenger mass can be 
inferred. Indeed the discovery of a higgs-hke boson in the mass region around 125-126 GeV by the CMS 
and ATLAS teams has important imphcations for supersymmetry. Notably, this mass is relatively large 
for the MSSM and large radiative corrections from stop/top loops are needed. Such contributions can arise 
through stops in the 0(10) TeV mass range or through lighter stops with maximal mixing [JHl ISHl HI]- 
Naturalness promotes the second option (recent discussions on the fine tuning issues can be found in 
ref. [22l [52]). A heavy Higgs with light stops can be obtained for large A-terms. In the absence of 
additional interactions ^-terms are generated through the renormalization group equation of the MSSM 
driven predominantly by the gluino mass. This requires a large messenger scale Mmess ^ 10"'^*' GeV and 
a heavy gluino rrig > 3 TeV [20j . Hence, a minimal model of gauge mediation seems to suggest heavy 
messengers. According to the figure 8 the gravitino for heavy messengers can be the dark matter of 
the universe if a late entropy production, either by the messengers or the spurion, has taken place. It is 
interesting to notice that for such large messenger scale the dilution caused by the spurion decay takes place 
automatically. Also, a heavy gluino does not modify the initial (before the dilution) gravitino abundance 
which is thermalized - given that Tih > Afmoss- 

On the other hand, it is possible to generate large A-terms at low messenger scales through superpo- 
tential interactions between MSSM and messenger superfields [211 US] . Such interactions generate sizable 
^-terms already at the messenger scale while not generating over-large one-loop rn?fj soft masses for the 
Higgses. This happens by taking into account a one-loop negative contribution to m\ that scales as 
~ A/Mtness- For low messenger scale, Mmess ^ A, this contribution is importrant. A direct cosmological 
implication of this MSSM-messenger mixing is that it renders the messengers rather shortlived (Ty ^ 0) 
to dominate the energy density of the universe. Moreover, the low messenger scale opens the parameter 
space towards smaller values of Afc where the thermal cross section of goldstino production from messenger 
thermal scatterings decreases (22). For sufficiently small coupling A and parameter k the Yj^'^^^'^' can 
be subdominant and, due to the MSSM-messenger mixing the goldstino yield from the messenger decays 
Y^,'^'^^ '^'^' is negligible. Given the cosmological constraint Vt^i2h^ < 0.11, the gravitino will be basically 
produced from scatterings with the MSSM plasma at the temperature Tf; where the i?-symmetry breaks 
-given that the helicity ±3/2 gravitino component is underpopulated, see figures 3 and 5. This is an ex- 
ample of gravitino dark matter production due to a symmetry violation. In addition, once messenger scale 
is indicated then the A and k parameters could be also probed. 

Other models exist that may lead to different conclusions. It is not the purpose of this work to present 
an exhaustive study all the GMSB models that aim to fit the LHC data. Here we briefly alluded to some 
general features of the GMSB phenomenology in the LHC era and exposed how the collider data can be 
directly connected, reconciled or contrasted with the gravitino cosmology. 

Finally, we note that, according to the findings of this work, cosmology implies a different kind of 
relation between the GMSB messenger scale and the gravitino mass. For a universe reheated only once 
(e.g. from the inflaton field) to high temperatures then the messenger scale is low, Mmoss < 10^ GeV and 
the gravitino mass in the range 1 GeV < m^/2 ^ 100 GeV. Otherwise, the gravitino is overabundnat, see 
figure 5. A heavy gravitino makes the NLSP a collider stable particle -it actually decays outside the solar 
system. In the case that the spurion field produces late entropy then the gravitino can be lighter. Of course 
all these constraints apply for Trh > Af,„ess. The reheating temperature of the universe is suggested to be 
that high, Tih > Af,„ess, by two reasons: the thermal selection of the supersymmetry breaking minimum 
[6] and the theory of thermal leptogenesis. The gravitino mass and the NLSP lifetime have additonal, 
important implications for cosmology. 

6.2 Further cosmological constraints 

Large scale structure 

Cosmological observations strongly indicate that the dark matter of the universe is cold [7] . The gravitino 
is a light particle and if too light it can behave like warm dark matter. The gravitinos in the early universe 
are relativistic and their momentum redshifts like p ex a~^ . However, non-thermally produced gravitinos 
by the decay of the spurion may be still relativistic at late times, for their production takes place at 
low temperatures T ~ McV - GeV. The free-streeming lenght of gravitinos has to be consistent with the 
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observational bound A/s < 0(100) kpc from the Lyman a forest data [9]. 

A„.100kpcy (^__j (^jjj^) (^^^j ^ (69) 

A supplementary bound comes from the reionization epoch 7 . These bounds constrain the gravitino mass 
from being too small and can be applied on our results. According to the figures 5, 6 and 7 there is a 
parameter region (upper right corner, i.e. where Afc -^ 1) where the gravitino is produced non-thermally 
from the spurion decay. The low messenger scale, M^ncss ^ 10^ GeV, and the relatively large values of 
the parameters A and k correspond to a light gravitino. Hence, this part of the parameter space can be 
excluded as a viable choice due to the observational bound (|69|) . 
BBN 

On the other hand, a heavy gravitino although compatible with the large scale structure constraints it 
makes the tnlsp large ([55)1 rendering the NLSP decay potentially dangerous for ordinary nucleosynthesis 
predictions. The BBN constrains the tnlsp to be up to ~ 10^ sec if the final-state particles are photons 
and this upper bound decreases significantly, nearly 10 orders of magnitude, for hadronic NLSP decays. 
A heavy gravitino 1 GeV < m3/2 ^ 100 GeV is predicted when there is no entropy production other than 
the very initial one; differently, the gravitino relic density exceeds the observational bound ^^3/2 ^ 0.11 
(see figure 5). When the gravitino mass exceeds a few GeV then cosmology requires either a comparatively 
heavy NLSP or/and an NLSP with small hadronic branching ratio or/and with small relic abundnace. In 
gauge mediation the NLSP is generically a neutralino or a stau and in some particular cases a sneutrino. 
The cosmology of a heavy stable gravitino, 1 GeV < 7713/2 ^ 100 GeV, can be viable and indicates towards 
specific gauge mediation supersymmetry breaking schemes. A recent discussion can be found in the ref. 

m- 

Leptogenesis 

The fact that the goldstino abundance is nearly independent of the reheating temperature with its yield 
being dominated either by the temperature T^ or M^^ss allows the thermal leptogenesis scenario to be 
realized without violating the bound ^13/2/7^ < 0.11. In the leptogenesis scenarios the baryon asymmetry, 
rjB — 6 ^ lO^^*', is generated in the out-of-equilibrium decay of right handed neutrinos and sneutrinos in 
the very early universe |53[ 119] . The right handed neutrinos can be generated by scatterings in the thermal 
bath for 

Trh > Ml - 10^ GeV , (70) 

where Mi the mass of the lightest of the heavy neutrinos. In this work we demonstrated that such high 
temperatures are allowed in the GMSB scenarios even in the absence of late entropy production given that 
the messenger scale is low enough, namely, Mmcss < 10^ GeV. If, on the other hand, the messenger scale 
is not that low then leptogenesis can still take place while satisfying the dark matter density bound due 
to dilution caused by the decay of the spurion field. The insensitivity of the helicity ±1/2 gravitino relic 
density to the reheating temperatures for T^h > Mmoss can help in achieving the correct dilution magnitude 
for the gravitino abundance without the washing-out of the baryon asymmetry. 

7 Conclusions 

The gravitino is a hypothetical particle predicted by a well motivated theory the supersymmetry which 
currently is being tested at the LHC. In the gauge mediation supersymmetry breaking (GMSB) theories the 
gravitino is the lightest supersymmetric particle and if i?-parity is conserved then it can be part of the dark 
matter of the universe. The gravitino relic density depends on its couplings and the reheating temperature 
of the universe. For reheating temperatures below the messenger scale the gravitino is produced solely by 
the MSSM particles. For reheating temperatures above the messenger scale the messengers have a dominant 
contribution and we actually find that the gravitino can easily attain a thermal equilibrium distribution. 
High reheating temperatures generically restore the C/(1)_r symmetry of GMSB sector suppressing the 
gravitino production from the MSSM superfields. Furthermore, the GMSB spurion X acquires a thermal 
vev about the origin of the field space and if its oscillations amplitude is large enough it can produce 
entropy at late times. 

The gravitino can have a relic abundance r23 /2 < JIdm for different values of the parameter space 
which controls the gravitino yield. If the messenger scale is low, lO^GeV < M,„oss < 10® GeV, and the 
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messenger fields have superpotential couplings with the MSSM superfields then the gravitino yield from the 
thermalized messenger scatterings and decays may not be able to saturate the 0dm bound. This is the case 
when GMSB messenger superpotential coupling A is sufficiently small. This is actually suggested by the 
conditions of the thermal selection of the supersymmetry breaking vacuum [5]. Moreover, the gravitino has 
to be heavy, 10 GeV < TO3/2 < 0(100) GeV hence, we conclude that there is an extra hidden sector, Z, that 
dominantly contributes to the supersymmetry breaking but has a subdominant, gravitational contribution 
to the soft masses. It is natural to expect that this hidden sector is characterized by an exact U(l)i{ 
symmetry |10j . The i?-symmetry, on the one hand, makes the .Z- vacuum an enhanced symmetry point 
and, on the other, implies that the Z-sector contributes mainly to the sfermion soft masses. 

The inadequate gravitino production from the thermalized messengers with Mmoss < 10^ GeV implies 
that the gravitinos are produced mainly by the MSSM sector. The MSSM superfields generate gravitinos 
dominantly at the temperature, T^, which is the characteristic temperature of the transition to the R- 
violating phase [18]. The gravitino relic density can have the observed dark matter density, 03/2/1^ — 0.11, 
in a quite constrained but not negligible part of the parameter space. For that low messenger scales the 
GMSB spurion can reach the zero temperature minimum without producing late entropy. 

When the messenger scale is large the gravitino yield from the thermalized messengers increases. The 
gravitino is thermally over-produced. However, a large messenger scale implies a large vev for the GMSB 
spurion, Xq, and it is found that generically the spurion field dilutes the thermal plasma and produces 
non-thermal gravitinos. It is possible the dilution and the non-thermal gravitino yield to have an acceptable 
size [31]. 

It is interesting that the current GMSB phenomenological models discuss the cases of low and large 
messenger scale. The input from particle phenomenology is essential for describing the cosmology of the 
gravitino in a more concrete and consistent way. In this work we have shown that gravitino cosmology can 
be compatible with high reheating temperatures with the upper bound imposed by the yield of the helicity 
±3/2 gravitino component. Therefore, the thermal selection of the supersymmetry breaking vacuum [6] 
and the leptogenesis scenario can be realized without facing the gravitino overproduction problem. The 
following years of the LHC data will hopefully shed "light" on the physics of dark matter. 
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A Thermal restoration of the U{\)^ symmetry 

A.l Exact i?-symmetry 

The connection between a global [/(l)i^ symmetry and the supersymmetry breaking is well established |10j : 
i?-symmetry for generic superpotentials is a necessary condition for supersymmetry breaking in the true 
vacuum and a spontaneously broken i?-symmetry a sufficient one. We expect i?-symmetry to be violated 
by different sources that may restore supersymmetry but not in a nearby vacuum. These sources may 
be 1/Mpi suppressed dimension- five operators in the superpotential or a constant term that cancels the 
cosmological constant in the vacuum. 

In the global supersymmetric limit the [/(!)/? symmetry can be an exact symmetry. A model that im- 
plements this behaviour has the minimal superpotntial (|T]) and Kalher the ([2]) with 64 = 1 and £5 < 0. This 
Kahler potential can originate from an i?-symmetric O'Raifeartaigh-like sector that breaks spontaneously 
the i?-symmetry, see for example ref. |14| . The _R-symmetry is a symmetry of the vacuum when X = 
and breaks spontaneously when X ^ 0. It is restored due to thermal effects at the temperature [5] [18] 
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For the minimal case of a 5 + 5 messenger sector, i.e. the </>, messenger quarks and leptons form a single 
complete SU{5) representation, it is A^ = 5. We can have additional SU{5) multiplets that couple to the 
spurion X field preserving the gauge unification. However, the gauge coupling at the GUT scale increases 
and if we require it to remain perturbative then we may add finite 5 + 5 pairs or a single 10 + 10 pair or 
(5 + 5)+(10 + 10) pairs to the particle content of the minimal SU{5) GUT. For large messenger scales 
values for the N as large as ~ 50 are allowed [8]. The Tfi can be quickly estimated: the spurion has 
negative squared mass 4F^/Aj at the origin and receives thermal corrections NX^T'^ /A from the messenger 
fields. 

A. 2 Approximate i?-syminetry 

Two basic examples of gauge mediation are the following [12l [13] : 

W ^FX + \X4)4> + c, W = FX + \X(t)4) - M(j)4> (72) 

with Kahler, K — \X\'^ — \X\'^ / Kl and N number of messenger fields (j) and in the fundamental represen- 
tation. Although these theories break the U{l)ii explicitly they are approximately _R-symmetric at high 
temperatures. Indeed, for high enough temperatures the i?-violating terms c and Mcjxj) are negligible and 
an approximate i?-symmetry restoration takes place. This can be seen from the evolution of the thermal 
average value for the i?-charged X field [5] : 

Tmn\ ) ^p2 I JV^2y2 ' min\ ) „p2 , jV^2y2 ^ ' 

The i?-symmctry breaking scale is the vev {X) = Xq. For the first case of gravitational stabilization 
the vev is the Xq = ch? /{2FMp^). For the second model, after the translation X — > M/\ — X it is 

Xq = M/\ thus, the Xq represents the scale of i?-breaking. 

According to (j73|) the thermal average value tends to restore the i?-symmetry. We can parametrize the 
degree of the i?-symmetry breaking by defining the parameter 6/j: 

6^(T)^ffil. (74) 

Temperatures higher than the cut-off scale are not expected (for A, > l{)^^Mp\) since a thermal equilibrium 
cannot be achieved. Thus, for the case of gravitational stabilization the second term at the numerator ()73|) 
is negligible. The parameter 5^ is given, for both cases, from the expression 

f4-^f 
bR{T) = \^ (75) 

Utilizing the definition ([7T|) of the Tfi we can recast the ([751) i^ito the simpler form 

hn{T) = ^--^. (76) 

Obviously, when T — )• the i?-symmetry breaking scale takes its maximum value, i.e. the zero temperature 
one, and when T — > oo the i?-symmetry is restored. In other words, the 6_r(T) parametrizes the i?-symmetry 
breaking scale at finite temperature with respect to the zero temperature scale. For the case of spontaneous 
breakdown of the i?-symmetry, discussed in the previously, the parameter bii(T) takes, approximately, the 
discrete values: 

bR{T>Tf^)^0 and ^^(r < T^) = 1 . (77) 

From the expression ([TG)) we see that bp){Tf)) ~ 0.5. Also, from ([TG)) we can re-derive the temperatures 
that the supersymmetry breaking vacua form, Tx , firstly given in ref. [5 , 6 . The Tx controls the thermal 
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selection of the supersymmetry breaking vacuum. It corresponds to the temperature that the minimum at 
the X-direction crosses the tachyonic boundary X = y^F/X. Hence, 



X{Tx)^ /f/A = bR{Tx)Xo (78) 

which gives the fohowing vahies for the parameter bj^; 

MT^)-2^^yy, hdTx )-^^- (79) 

It has to be bB.{T^x^), hniT'^x^^) < 1 and from ^ we take 

9 8 c /f , ,16 FM [f , , 

Tl-^ ^\ — and T|- ~ tt^a/— (80) 

^ N AAf 2j V A ^ N X^Al V A ^ ' 

which are the temperatures derived in ref . [S] [6] . We also note that Tx > Tfi . 

B Gravitino production from the messenger sector 

The gravitino yield ^3/2, i-e. the number to entropy ratio, is found by solving the Boltzmann equation 
describing the gravitino production 



dT sHT 



\ i ij J 

The scattering processes of messenger fields contributing to goldstino production are X^ + ^^ ^^ </" + V'G ; 
4> + X^^V^ + ^G, X0 + A° ^ X0 + ^G, (/) + A" ^ </.* + Vg, X0 + X0 ^ A'^ + ^G and + </>* ^ A« + ?/;g. 
The corresponding cross sections to goldstinos given in ref. [5S] read 



^ ,, ., , ^ ^ ., 2(2,s2-3sAf2_ + M4^^j + ,(5_2A/^_)log(^) 

2^ cr(A,„css + BmSSM -^ Ancss + V'g) = CA fc ( _ ]\f2 \2 ^^^^^^ 

A.B,B' '^^'^ ^^VossJ 

(82) 
and 

1 V^ / . ./ r. , X ^x2,2 ^^ "t^^^mcssJJ + JW^egglOg I 2 j^ 2 .ji/2 

T: E ^(^mess+4ncss ^ ^^MSSM+^g) ^ CA'^^" .^_.,.2 A ""^""^^ 

(83) 
where ^moss — <f>, X<t> and Bmssm = A", V^ are respectively messenger and MSSM gauge superfield com- 
ponents. 

The decay width of a spin 1/2 messenger to a scalar messenger plus a goldstino due to the Lagrangian 
interaction 5C = Xk ^g Xf 4' reads [3S] 

p^2 / M^V 

r.(x,- ^ 0^g) = 4 X -^M,^ y-i4:) ■ ^''^ 

The fermion messenger mass is Af^- — XXq = Afmcss and the squared mass of the scalar messengers is 
M| = M^ggg ± XF. Hence, the decay rate is recast into 

r.(x^^0^G) = ^AV^. (85) 

mess 

Messenger Scatterings 

The thermally averaged product appearing in (j8ip is found in 35 and takes the form 

(c^«+i^V'G+... % '^i '^i) = FT / cis Vs ^1 ( ^TT ) Ipl^crtot(s) (86) 
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Figure 8: The scattering cross sections (j82|) and (j83| of thermalized messengers to goldstinos as a function of the 
center of mass energy squared. The figures show the fast growth of the cross section close to the messenger mass 
squared and the 1/s fall at high energies. Here the exemplified messenger mass is Af,„css = 10^ GeV. 



with (Ttot the sum of the scattering cross sections (1821) and ((83|) . Ki the modified Bessel function of the 
second kind of degree one, and \p\ the absolute value of the incident i particle in the center of mass frame. 
The scattering yield is given by the expression 



Y, 



mess(sc) 
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with X — {Mi + Mj)/T. When one of the incident particles is a messenger then Mi — M^icss and M2 — 
we have the yield 
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Gravitinos reach equilibrium abundnace when Tg^ > H, see eq. (pij) . for at least one expansion time. The 
relic abundance of thermal gravitinos is given by the expression |35j 



3/2, 



0.278 X 3/2 



1 — exp 



^'*-dT' 



HT 



(89) 



Messenger Decays 

The gravitino yield from the decays of the thermalized messengers is given by 
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where the sum is over all the messenger particles. The yield is dominated by temperatures T ^ M„ 
The decay width of messengers to goldstinos is given by the expression 



rA(x0^0V'G) = 4x x^k 
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When it is T^'Jf" = Tx the yield ^ takes the value 
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(92) 



30 



where Nxi the number of the messengers that decay to hghter messenger superpartners plus the goldstino. 
It corresponds to the rehc abundance 

However, this goldstino yield from messengers decays changes when there are more decay channels. Actu- 
ally, the goldstino yield from messenger decays cannot exceed the thermal equilibrium value of relativistic 
messenger fields, Y'^^N^^ ^ 1Q~^N^^, because the goldstino is one of the two daughter particles from 
the decay of the messenger, the mother particle. Moreover, given the existence of other decay channels 
of messengers to MSSM particles with decay width Ty and r„j the goldstino yield from messenger sector 
decays is bound to be less than the thermal equilibrium value regardless the value of the product A • k 
and the messenger mass M^css- Taking into account the other decay channels of the messenger fields the 
goldstino yield from the messenger decays is 

Y--^^^^^^B,„N,^Y^^^^ (94) 

where B3/2 = rA/(rj, + T™ + Ta + ...). Furthermore it is K{T = M,„ess) = ^tot/H » 1. 

The expression (^^ has a Af "^^^ dependence and appears to dominate over the goldstino yield from 
scatterings, ([55)1 . for Mmess ^ lO^GeV. However, when there are other messenger decay channels this con- 
clusion is altered. In particular for Afc < 10^^ the product Afc, according to the ratios (02]) and (03]) depicted 
in the figure (14.2.21) . render the decay rate T\ negligible compared to the Ty and Tm- Also, the goldstinos 
do have a thermal equilibrium abundance from messenger scatterings for Afc > 10~^(Mmoss/10*GcV)^/2 
(see relation ([32])). Hence, we find that B3/2 ~ rA/(rj^ + T„,) < {10^^^ - 10"^} and it is 

T^mcss(dcc) T^mcss(sc) ^ -^eq /q-n 

^3/2 "^ ^3/2 - -^3/2 ' y^"^) 

for any messenger scale. 
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